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RANDOM POLYNOMIALS WITH PRESCRIBED NEWTON POLYTOPE 


BERNARD SHIFFMAN AND STEVE ZELDITCH 


Abstract. The Newton polytope Pf of a polynomial / is well known to have a strong impact on its 
behavior. The Bernstein-Kouchnirenko theorem asserts that even the number of simultaneous zeros in (C*)^ 
of a system of m polynomials depends on their Newton poly topes. In this article, we show that Newton 
polytopes further have a strong impact on the distribution of zeros and pointwise norms of polynomials, 
the basic theme being that Newton polytopes determine allowed and forbidden regions in (C*)"^ for these 
distributions. 

Our results are statistical and asymptotic in the degree of the polynomials. We equip the space of 
polynomials of degree < p in m complex variables with its usual SU(m + l)-invariant Gaussian probability 
measure and then consider the conditional measure induced on the subspace of polynomials with fixed 
Newton polytope P. We then determine the asymptotics of the conditional expectation of 

simultaneous zeros of k polynomials with Newton polytope NP as A —> oo. When P = S, the unit simplex, 
it is clear that the expected zero distributions E|;vi:(^/i,...,/fe) uniform relative to the Fubini-Study form. 
For a convex polytope P C pS, we show that there is an allowed region on which is 

asymptotically uniform as the scaling factor N oo. However, the zeros have an exotic distribution in the 
complementary forbidden region and when k = m (the case of the Bernstein-Kouchnirenko theorem), the 
expected percentage of simultaneous zeros in the forbidden region approaches 0 as A ^ oo. 
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1. Introduction 


It is well known that the Newton polytope P of a polynomial f{zi,...,Zm) of degree p has a crucial 
influence on its value distribution and in particular on its zero set. Even the number of simultaneous zeros in 
(C*)™ := (C\ {O})’" of m generic polynomials fi,..., fm depends on their Newton polytopes |Be[IKolllKo2| . 
Our purpose in this paper is to demonstrate that the Newton polytope of a polynomial / also has a crucial 
influence on its mass density \f{z)\'^dV and on the spatial distribution of zeros {/ = 0}. We will show that 
there is a classically allowed region where the mass almost surely concentrates and a classically forbidden 
region where it almost surely is exponentially decaying. The classically allowed region is the inverse image 
of the (scaled, open) polytope ^P° under the standard moment map p of complex projective space 
CP™. The simultaneous zeros in (C*)™ of m generic polynomials /i,..., /„ with Newton polytope P tend 
to concentrate (in the limit of high degrees) in the classically allowed region, giving a kind of quantitative 
localized version of the Bernstein-Kouchnirenko theorem 1^ 1^ IT^ . Polynomials with a given Newton 
polytope P are often called sparse in the literature, and methods of algebraic (including toric) geometry have 
recently been applied to the computational problem of locating zeros of systems of such sparse polynomials 
(e.g., see EHl Ell EH E31^). Our results give information on the expected location of zeros when the 
polynomials are given the conditional measure described below. To our knowledge, this kind of precise 
asymptotic concentration of zeros of sparse systems of random polynomials has not been observed before. 

The Newton polytope has an equally strong (though different) impact on the common zeros oi k < m 
polynomials /i,..., /^. The image of the zero set of /i,..., /^ under the moment map is (up to a logarithmic 
re-parameterization) known as an amoeba in the sense of [HkzIImEI- Results on the expected distribution 
of amoebas can be obtained from our results on the expected zero current; for example, for a polytope ^P 
with vertices in the interior of the standard unit simplex E C K™, there is also a very forbidden region which 
the amoeba almost surely avoids (see Corollary II.till . 


1.1. Statement of results. The patterns in zeros discussed here are statistical—they hold for random 
polynomials with prescribed Newton polytope—and are asymptotic as the degree of the polynomials tends 
to infinity. To state our problems and results precisely, let us recall some definitions. Let 

/(zi, . . . , Zm) = 'y ' aaXct (-21, . ■ . , Zm}i Xa (z) ~ ' ' ' ^rn' (^) 

(a = (ai,..., am), |<t| = ai -|- • —h am) be a polynomial of degree p in m complex variables. By the support 
of / we mean the set 

5/ = {a G Z™ : aa 7 ^ 0}, (2) 

and by its Newton polytope Pf we mean the set 

Pf := the convex hull in R.™ of Sj. (3) 

Our aim is to study the statistical patterns in the space 

Poly(P) = {/ G C[0i, ...,Zm]:PfCP} (4) 


of polynomials with support contained in a fixed Newton polytope P. We therefore need to find a natural 
measure on this space of polynomials. Since our purpose is to compare zero sets and masses as the polytope P 
varies, we view the polytope P as placing a condition on the Gaussian ensemble Poly(pE) of all polynomials 
of degree < p (where pS is the dilation by p of the standard unit simplex S), and we give Poly(P) the 
resulting conditional probability measure. Here, p is chosen such that pE ID P. 

To describe this Gaussian ensemble, we first identify Poly(pE) with the space of homogeneous polynomials 
of degree p in m-|- 1 variables in the usual way, i.e., by identifying / G Poly(pE) with the homogeneous poly¬ 
nomial F such that P(l, zi,. ■ ■, Zm) = f{zi ,..., Zm)- Using this identification, we give the space Poly(pE) 
the SU(m -(- l)-invariant inner product 


(/l,/2) 


(Pi, P2 )s2"*+i 


1 

to ! 


/S2m+1 


P 1 P 2 diz , 


( 5 ) 


where dv is Haar probability measure on the {2m + l)-sphere C 0™+^. The monomials Xa are 

orthogonal, so an orthonormal basis of Poly(pE) is given by {||Xa||~^Xa}|„|<p, where || • || denotes the norm 
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in Poly(pE) given by lO- We emphasize here that the norms {Hxali} of the monomials {xa} do not depend 
on the constraining polytope P —they are given by the SU(to + l)-invariant inner product ISJ on Poly(pI]). 
(Of course they do depend on the choice of the integer p.) 

The SU(77i + l)-invariant Gaussian measure Xp corresponding to the inner product @ is defined by 


dlpif) 



/=E 

|a|<p 



( 6 ) 


where kp = #{a S N™ : |a| < p} = {"^p^)- Thus, the coefficients Cq are independent complex Gaussian 
random variables with mean zero and variance one. We then endow the space Poly(P) with the associated 
conditional probability measure Xp\P'- 


dlp\p{f) 


-^e Id'dc, / = ^ ( 

aeP 


Xo 


llxo 


(7) 


where the coefficients Cq are again independent complex Gaussian random variables with mean zero and 
variance one. (By an abuse of notation, we let J2aeP denote the sum over the lattice points a G P H Z™; 

denotes the cardinality of P n Z"*.) We observe that Poly(P) inherits the inner product (/i,/ 2 ) from 
Poly(pI]), and that 7 |p is the induced Gaussian measure. Probabilities relative to yjp can be considered as 
conditional probabilities; i.e. for any event P, 


Prob.^{/ e E\Pf =P} = Prob.^|p(P). 


The expected distribution of mass and zeros of polynomials with fixed Newton polytope turns out to 
involve the moment map p, : (C*)™ ^ K™ given by 


p{z) = 


l^il 


,1 + lkP’ ’1 + lkP 

The image of p is the interior of the standard unit simplex S in MP with vertices at the points 


( 8 ) 


( 0 ,..., 0 ), ( 1 , 0 ,..., 0 ), ( 0 , 1 ,..., 0 ), ...,( 0 ,..., 0 , 1 ) . 


The map p is the moment map of CP™ D (C*)™ (with the Fubini-Study symplectic form ccps) and plays a 
role in the geometric approach in IS^ IST^ f STp . where we regard polynomials with a fixed polytope as 
sections of a holomorphic line bundle on a toric variety (see 

By a convex integral polytope P, we mean the convex hull in R™ of a finite set of points in Z™. We use 
the moment map p to describe our classically allowed regions: 


Definition: Let P be a convex integral polytope in K™ such that P C pS. The classically allowed region 
for polynomials in Poly(P) is the set 

Ap ■-p-^ (^P°^ C (9) 

(where P° denotes the interior of P), and the classically forbidden region is its complement (C*)™ \ylp. (If 
P has empty interior, i.e. if dimP < m, then Ap = 0.) 


All of our results on the asymptotic expected distribution of zeros of systems of random polynomials with 
prescribed Newton polytopes follow from our asymptotic result ('Theorem 11.211 on the pointwise expected 
values of random polynomials with given Newton polytopes. Our results on zeros have a somewhat different 
flavor as the codimension (or number of polynomials in the system) varies. We first describe the results for the 
maximum codimension m where the zero set is almost surely discrete ITheorem ll.lll . Then after describing 
the expected value distribution, we discuss the codimension 1 case, where the zero set is a hypersurface, and 
we finish with a result for general codimension that encompasses the previous results. 
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1.1.1. Distribution of zeros: the point case. Let us first consider the simultaneous zero set of m independent 
random polynomials in m variables. Bezout’s theorem tells us that m generic homogeneous polynomials 
Fi,..., Fm of degree p have exactly p"* simultaneous zeros in (C*)"*. In fact, one immediately sees (by 
uniqueness of Haar probability measure on CP™) that the expected distribution of zeros is uniform with 
respect to the SU(TO+l)-invariant projective volume form, when the ensemble is given the SU(m+l)-invariant 
measure d'jp. 

According to the Bernstein-Kouchnirenko theorem the number of common zeros in (C*)™ 

of m generic polynomials {/i,..-,/m} with given Newton polytope P equals m!Vol(P), where Vol(P) is 
the Euclidean volume of P. For example, if P = pE, where E is the standard unit simplex in R™, then 
Vol(pE) = p™Vol(E) = and we get Bezout’s theorem. (More generally, the fj may have different Newton 
polytopes, in which case the number of zeros is given by the Bernstein-Kouchnirenko formula as a ‘mixed 
volume ’ |BelIKollIKo2| .i 

Now consider the ensemble of m independent random polynomials with Newton polytope P, equipped 
with the conditional probability 0 . We let E|p(Z/^^,,. denote the expected density of their simultaneous 
zeros. More generally, for TV S we define the expected zero density measure E| 7 vp(^/i,.by 

^\Np{Zfi,...J^){U) = J d7jVp|Arp(/l) ■ ■ ■ y d-'yNp\Npifm) \ff{z € C/ : fi{z) = ■ ■ ■ = fmiz) = 0}] , (10) 

for U C (C*)™, where the integrals are over Poly(fVP). In fact, E|^p(Zjj^...j^) is an absolutely continuous 
measure given by a C°° density (see Proposition 15.111 . 

Our first result shows that, as the polytope P expands, these zeros are expected to concentrate in the 
classically allowed region and have (asymptotically) uniform density there. We measure the density of zeros 
with respect to the projective volume dVolcp™ = where tups = ^^'91og(l + ||^P) is the Fubini-Study 

Kahler form on C™ C CP™. (Note that with this normalization, the volume of CP™ is ^.) 


Theorem 1 . 1 . Suppose that P C pE C K™ is a convex integral polytope with nonempty interior. Then 


lim 
AT—^oo 


1 

(iVp)™ 


E|Afp(^/l...../m) 


iv™g on Ap 

0 on (C*)™\Ap 


in the measure sense; i.e., for any Borel set B C (C*)™, we have 

^^^^E|Arp (#{0 G B : fi{z) = ■■■ = fm{z) = 0}) ^ mlVolcp™ (P H Ap) . 

Theorem ll.ll is a special case of our general result on zeros fTheorem ll.411 . In fact, our results imply that 
the convergence of the zero current on the classically allowed region is exponentially fast in the sense that 

^\Np{Zf,^...jJ = (iVp)™u;^s + O (e-^^) on Ap , 
for some positive continuous function A on Ap. 

Note that we can choose any integer p such that pE D P. Theorem 1 (as well as our other results) depends 
quite strongly on the choice of p, since the classically allowed region Ap shrinks as p is increased. In order 
to minimize notational clutter, p does not appear Ap, E|jvp, |/(2 )|fs, and some other expressions to be 
define below. 

Figure ^ shows the classically allowed region (shaded) and the classically forbidden region (unshaded) 
when P is the unit square in (and p = 2): 


1.1.2. Mass distribution. By the mass density of a polynomial / G Poly(pE) at a point z € (C*)™, we mean 
the square |/( 2 )|pg of the Fubini-Study norm 


I/WIfs := 


lf(z)l 


(1 


2)p/2 


= c = 




= (l,Zi,...,Z,„)ep2m+l 


( 11 ) 


where F is the homogenization of / described above. Our next result describes the asymptotics of the 
expected mass density of £^-normalized polynomials with Newton polytope NP, i.e. the expected den¬ 
sity Ey„p(|/(z)||g) with respect to the Haar probability measure, denoted by p/vp, on the unit sphere in 
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^2! 




Figure 1. The classically allowed region for P = [0,1] x [0,1] 


Poly(-/VP). We shall show that Ej,^p(|/( 2 ;)||g) is asymptotically uniform with respect to Fubini-Study mea¬ 
sure in the classically allowed region (as if there were no constraint at all); while in the forbidden region 
the mass decays exponentially. Thus in the semiclassical limit ^ 00 , all the mass concentrates in the 
classically allowed region. 

Behavior of the expected mass density (as well as the zero distribution in positive dimension) in the 
forbidden region is subtle, so we pause to introduce the relevant concepts. We show in non that the 
classically forbidden region decomposes into the disjoint union of the normal ‘flow-outs’ of boundary points 
of P, 

Flow(a;) := {e^ ■ z : t G Cx, PfJ-iz) = x) , x G dP n pE° , (12) 

where Cx C M"* is the normal cone to P at a; (see il2.1ll . Here r ■ z = {nzi,, rmZm) denotes the RIJ* action 
on (C*)™. 

We use these flow-outs to divide the forbidden region into subregions TZp as follows: Recall that the 
boundary of a polytope P decomposes into a disjoint union of its faces (of all dimensions). For each face F 
of P not contained in 9(pE), we let 

TZf = [J Flow(a:) = {e’’ ■ z : t G Cp, pp{z) G P} . (13) 

xeF 

where Cp is the normal cone to P along the face P; i.e., the normal cone to P at any point of P (see O- 
Note that for the case where P is the open face P°, TZp is the classically allowed region Ap. We call TZp the 
flow-out of the face P. The regions TZp all have nonempty interior. A sample illustration of these regions if 
given in Figure 01 in 

We can now state our result on mass asymptotics: 

Theorem 1.2. Suppose that P C pT, C K.™ is a convex integral polytope such that P 5(pS). Then 
the expected mass density of random -normalized polynomials with Newton polytope NP is given by the 
asymptotic formulas: 

Ei/„p (|/(^)||s) ~ Co-\-ciN~^-\-C 2 N~^-\ -, for z G Ap, 

(|/(z)||s) = [c^(^) + o(iv-i)] , for z G TZf, F C , 

where cq = yo\(p) ’ ^ ~ codimP, Cq G C°°{TZf), and bp is a positive function on {C*)'^\Ap. Furthermore, 
the remainder estimates are uniform on compact subsets of the open regions TZf and of Ap. 
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We shall give a formula for bp{z) below. A more precise asymptotic formula for the expected mass density 
is given by Theorem 14.II isee also (O). 

It follows that 


E.„p (|/(^)||s) 


Vol(P) ’ 

0 , 


for z G Ap 
forzG (C*)’"\:4^ 


as illustrated in Figure |21 below (plotted using Maple) for the case where P is the unit square. (Recall Figure 
^for the depiction of Ap for this case.) 


N=10 


N=100 




Figure 2. (l/(^)lls) ^r P = [0,1] x [0,1] 


Heuristically, the mass concentration in Ap can be understood as follows: the mass of a monomial Xa G 
Poly(A'p) concentrates (exponentially) on the torus The constraint Sf C NP thus concentrates 

all the mass in Ap, with exponentially small errors, and the mass there is uniformly distributed with small 
errors. Such mass concentration of monomials is given in formula and in more detail in (Sill]. 

We also show that logE|jvp (|/(2 :)Ifs) ^ —bp{z) uniformly on compact subsets of (C*)™ (Propo¬ 
sition Oil . even at the interfaces of the regions TZp. Combined with the Poincare-Lelong formula, this 
uniform convergence result leads to the asymptotics of the expected distribution of zeros of a single random 
polynomial with polytope NP. Using independence of sections and the Bedford-Taylor Theorem Ell on 
continuity of the Monge-Ampere operator, we obtain asymptotics for any number 1 < A: < to of sections. 

We shall show that bp is C°° on the regions TZp, but is not even on the interfaces between these regions. 
To state our formula for bp, we let q(z) denote the unique point in dP such that z G Flow( 9 (z)) and we let 
Tz G Cq(z-) be such that z = ^ where ^ G ^ig(z)^. Thus g(z), are given by the conditions: 

• z) = q{z) G dP, (14) 

Tz e (15) 

The existence and uniqueness of q{z) and Tz are stated in Lemma 14.31 
We then have: 

f 1 + lkf \ 

\i+\\e--^/^-zrj 


bpiz) = -{q{z),Tz) +p\og 


for zG (e)™\Ap . 


(16) 













RANDOM POLYNOMIALS WITH PRESCRIBED NEWTON POLYTOPE 


7 


To understand formula better, for any x G P (not necessarily a lattice point), we consider the 
‘monomials’ 


|X,(Z)|:= |zr = |zi| 


> \Xx{z)\ ■■= |Xa:(7:)|FS = 


\Xx{z)\ 


(l+||zP)P/2 


7 .- 


\Xx{z)\ 

supixxl 


(17) 


for z G (C*)™. Thus the normalized monomial Mx has sup-norm 1; in fact it takes its maximum on the 
torus ^“^(ix). One easily checks that 

bp{z) = -21og7W,(^)(z) , (18) 

and thus Theorem o says that 

E..P (l/(^)lls) = [co (^) + 0{N-^)] N-^/^ M,(,)(zf^ , zG 7Z°p. (19) 

We will also obtain an integral formula for bp (see Proposition l5.dll . which allows us to interpret bp as an 
‘Agmon distance’ (see HI.411 . and q{z) is the closest point of P to pfi{z) in this sense. 

We do not assume in Theorem o that P has interior. However, if P° = 0 (i.e., if dimP < m), then we 
must assume that P ^ d{pYj) since if P is contained in a face of pS, the decay function bp is not defined. 
Note that if P° = 0, then Ap is vacuous and the first expansion of the theorem does not occur. 


1.1.3. Distribution of zeros: arbitrary codimension. Our next results concern the zero set of one or more 
random polynomials. Since the zero set 

\Zfu...,h h= e (C*)- : /i(z) = ... = /fe(z) = 0} (20) 

is a submanifold of complex codimension k (without multiplicity), for almost all polynomials /i,..., /fc, it 
defines a current of integration G We recall that this current is given by 

[ p, for test forms p e ^"^-'‘’■'"-^(C*)™) . (21) 

PZh ./J 

The zero set P/i,...,/fc| also carries a natural (2m — 2fc)-dimensional Riemannian volume measure (induced 
from the Fubini-Study metric on CP"*), denoted ||, against which one can integrate scalar functions. 

Since the volume form on any complex u-dimensional holomorphic submanifold of CP"* is given by the 
restriction of ^Wpg, the volume measure on the zero set is given by = ^/i,..../fc A 

i.e., 

iWZfu-JkW,^) = [ (pdVohm- 2 k = T’Wps"'', for p e X>((C*)"*) . (22) 

As a measure, ||^/i,..../fc||(H) = Vol(|Z/^^,,.| n U) for open sets U C (C*)"*. (For the theory of currents 
defined by complex algebraic or analytic varieties, see for example, |Shl 1.3].) We shall discuss the expected 
values of both the current of integration . and the measure 

We first consider the expected distribution of zeros of 1 polynomial. We denote by E|p(Zj) = 'E.y^\p{Zf) 
the conditional expectation of the zero current of a random polynomial / G Poly(P) with Newton polytope 
P. In fact, E|p(Z/) is actually a smooth (1, l)-form on (C*)"* (Proposition 

Let us recall what happens when P = pH. By the uniqueness of the SU(m-|- l)-invariant Kahler form wps 
on CP"*, the expected zero current 'E(Zf) taken over all polynomials of degree p is given by pwpSj where 
wps = ^ddlogWzW^ is the Fubini-Study Kahler form on CP"*. Thus the expected distribution of zeros, 
as well as the tangents to the zero varieties, is uniform over CP"*. We now describe how the expectation 
changes if we add the condition that Pf = P. 

Theorem 1.3. Let P C pH C K*" be a convex integral polytope such that P d(j)H). Then there exists a 
closed semipositive (1, l)-form ipp on (C*)"* with piecewise C°° coefficients such that: 

i) iV-iE|wp(Zy) ^ i/;p inCl^((C*r). 

ii) ipp = pujps on the classically allowed region p~^(^P°). 

iii) On each region PPp, the (1, l)-form ipp is C°° and has constant rank equal to dimE; in particular, 
if V G pH° is a vertex of P, then 'f’p\n° = 0. 
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We see from part (iii) that the zero set \Zf\ of a polynomial with polytope NP typically intersects the 
classically forbidden region \ ^P) in the semiclassical limit N —>■ oo. However, there are subtler ‘very 

forbidden regions’ that \Zf\ avoids in the case where the polytope has vertices in the interior ofpE, namely 
the regions TZ° comprising the flow-out of these vertices. 

As a corollary, we obtain some statistical results on the so-called ‘tentacles’ of amoebas in dimension 2 
(see ^<5.511 . Roughly speaking, the (compact) amoeba of a polynomial f{zi,Z 2 ) is the image of the Riemann 
surface Zf under the moment map ^ on (C*)^, and the tentacles are the ends of the amoeba. Certain 
tentacles must end at vertices of the triangle E while others are ‘free’ to end anywhere along the boundary 
of E. In Corollary 15.61 we will prove that (in the limit N oo) almost all of the free tentacles of typical 
amoebas tend to end in the classically allowed portion of 9E. 

We call the form ipp in Theorem 11 .di the limit expected zero current. Our explicit formula for -ipp is 

tpP = pwFS - ^^^^ddhp , (23) 

ZTT 

where bp is given by m- In fact, (|13 holds as an equation of currents, and the current tpp G I?'^’^((C*)™) 
is closed and positive. By convergence in (i), we mean P convergence of the coefficients on every 
compact subset of (C*)™. (Recall that E|jvp(^/i,...,/fc) is a {k, A:)-form with smooth coefficients.) If we write 
ipp = 1S V’jfc {z)dzj A dzk, then (jpjk^z)') is a semi-positive Hermitian matrix. By the rank of ipp at 

z, we mean the rank of the matrix (p:ji^{z)'). Note that if TZp and TZpi are adjoining regions (i.e., have a 
common codimension I interface), then F and F' are of different dimensions, so ipp must be discontinuous 
along the interface. 

Boundary points of the TZp are called transition points. The set of transition points comprises the 
discontinuities of ipp. Points of dAp are always transition points, and there will be others whenever P has 
a face F C pE° of codimension at least 2. For example, consider the case where P is the unit square. In this 
case, there are two interior faces whose flow-outs are the connected components of the classically forbidden 
region (see Figure^] above) and the set of transition points equals dAp. We shall also give an example with 
an interior vertex (see where the forbidden region is connected but decomposes into flow-outs of 3 

faces (see Figure 0)1 and there are transition points not in dAp. 

The form ipp not only encodes the expected (normalized) density of the zero set, but also the expected 
density of tangent directions to the zero set. In the course of the proof of Theorem II .31 we will show that 
in the forbidden region, the limit tangent directions are restricted. In particular, as the polytope expands, 
the tangent spaces to typical zero sets approach tangency to the ‘normal flow’ {e'’’’*'*® • : r, 6* G Tp C K"*}. 

A precise formulations of this fact is given in Theorem EB Thus, while the expected distribution of zero 
densities is absolutely continuous, the expected distribution of zero tangents is singular. 

Finally, we consider the general case of k < m independently chosen random polynomials 

fj G Poly(Pj) , Pj C PjE , Pj ^ 5(pjE) (1 < j < fc) , 

and we let Ejp^ p^(Z/j_...jj,) denote the expected zero current with respect to the probability measure 
7pi|Pi X • • • X 7pfc|pj, on the product space. If Pi = • ■ ■ = Pk = P, then we also write E|p p(Z/j_,..j^) = 

E|p(^/i,..../fc)- 

Theorem 1.4. Let Pi,..., Pk be convex integral polytopes in R>g. Then 

N-%^P,_^P,{Zf,_j,)^iPpp---AiPp, tnCUiCnn, as N^oo. 

Theorem II. II is a consequence of Theorems 1 1.31 and II. 41 Furthermore, the expected volume of the zero set 
of a system of k polynomials has the following exotic distribution law, as given by the following corollaries: 

Corollary 1.5. Let Pi,...,Pk be convex integral polytopes in R.>q. Then for every relatively compact, 
open set U C (C*)"*, we have 

^E|Arp^_..._jvpj^Vol(|%,...,/J n U) _ fc)! ^i’Pk A ■ 

Corollary 1.6. The expected zero current A^“^E| 7 vp(^/i,...,/fc) tends to 0 at all points of each forbidden 
subregion TZp with dimE < k. 
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1.2. Examples. We illustrate the results of Theorems in some simple cases in dimension m = 2. 

For our examples, we shall describe the forbidden subregions TZf given by itnii- fiH) and compute the decay 
function bp(z) and the limit expected zero current tpp. 

To simplify our computation of ijjp from (EHll , we write 

u{z) = -bp{z)+plog{l + ||zf) =plog(l + . ^|| 2 ) ^ , 

so that 


V'P = ^C^ddu = 


27T 


27T 


plog(l + ||e ■ zW^) + {q{z),T^) 


(24) 


1.2.1. Example 1: the square. For our first example, we let P be the unit square with vertices {(0,0), (1,0), (0,1), (1,1)} 
and we let p = 2. Recalling that 


1^{Z1,Z2) = 


Ri| 


F2| 


,l+klP + |z2P ’ 1+|Z1P + |Z2|V ’ 
we see that the classically allowed region is given by 

Ap = {izi,Z2) : |zip - 1 < \z2\‘^ < Izip + 1} , 
as illustrated in Figure ^ The forbidden region consists of two subregions: 

Tip = {izi,Z2) : |z2p > \zi\^ + 1} , F = {{xi, 1) : 0 < ti < 1} , 

Up* = {{zi,Z 2 ) ■■ |z 2 p < kip - 1} , F* = {( 1 , 2 : 2 ) : 0 < 2:2 < 1} ■ 

Suppose that z is a point in the upper forbidden region TZp. Recalling we write Tj = (ti , r 2 ); 

then Ti = 0 since r T Tp. The boundary point q{z) of P whose normal flow contains z is given by 

qiz) = 2 p{e-^F 2 . z) = {a,l) G F C dP . 


Writing 
we have 

Therefore 


|zip = Si, |z2p = S2 


Si 


|e ^^/^Z 2 p = e '^^52 = S 2 


S2 


1 + Si + S2 


Si = 


1 + Si + S2 


Si 


1 — a’ 


We have 


log(l 


^-rj2 . ||2n ^ 



1 

Si 

S2 


= —, a = 


1 — a 

a 

3—J-2 

= S 2 /S 2 ^ 

Si 1H 



as2 \ 

log ( 


: , 1 + kiP 
1^12 


Zl 


1 + Si l+|zip’ 
_ 1... |2 


^2 I 


{q(z),T^} = 


1^11 


1 + \zi 


| 2 > ■ 


1 , 0, log 


1^2 


^2 I 

2 


Z 2 p = log(l + Izip) +log2 


1 + Uil 


= log|z 2 p - log(l + IziH . 


Therefore 


= log|z 2 p +log(l + |zip) + log4 . 


We conclude that 

^aaiog(l + |zi|2) 


ipP = < 


(klP + l<k2n 


for z e TZf 

^^5i91og(l + |zip + |z 2 p) = 2a;FS for z G .Ap (|zip - 1 < |z 2 p < |zip + 1) 


. ^a51og(l + |z 2 p) for zGP^. (|z 2 p < |zip-l) 

(where the third case is by symmetry). Note that ipp has constant rank 1 in both of the forbidden regions 
TZp, TZp-, as indicated in Theorem a iii). 
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On TZf, 'vre have: 


4|z2p(l + kiP) 


for 


I-Zll^ + 1 < |Z2|^ 


(1+|ziP + |z2|2)2 

On the boundary {|zip = |z 2 p — 1}, we have = 1 as expected. On {\zi\ = c}, we have the growth 

rate ^ 1/|-Z2p as Z 2 oo. To obtain on we interchange zi and Z 2 in the above. 


1.2.2. Example 2: a trapezoid. We next consider the trapezoidal polytope of Figure |21 below. Comparing 
with the case of the square, we see that the classically allowed region is given by \z 2 \^ < + 1, and the 

forbidden region coincides with the upper forbidden region TZp = {|z 2 p > l^iP + 1} from Example 1. Thus, 
the map z (tz, q{z)) is the same as before when z is in the forbidden region TZp. Hence, and ipp 

are also the same as in Example 1 on TZp. On the classically allowed region, _ i _ 2a;ps. 



Eigure 3. P^I? = {(0,0), (1,0), (2,0), (0,1), (1,1)} 


1.2.3. Example 3: trapezoids of higher degree. Now let n > 2 and let P be the trapezoid with vertices 
(0,0), (n + 1, 0), (0,1), (1,1) given in Figure 0] below. In this case, p = n + 1 and P has an interior vertex 
V = (1,1) e s°. 

We see that the classically allowed region p~^{^P°) is given by 

= |(^i,22) : k2p < 

This time, the forbidden region consists of three subregions: Pp, Pv, Pp', where 


F = {{xi, 1) : 0 < xi < 1} , F' = {{xi,X 2 ) ■ X 2 = ^{n + 1 - xi), 1 < xi < n + 1} . 
(See Figure 0 below.) 

Suppose that z is a point in the region Pp. Then Tz = (0,r2), and 

q{z) = {n + l)/r(e“'^*/^ • z) = (a, 1) S (n + 1)F C dP (0 < a < 1) . 


Again writing 
we have 

Therefore 


|ziP=Sl, |Z2|^=S2 


|e-"^/2z2p=e-"=S2 = S2, 


Si 


S2 


1 


1 + Si + S2 71+1 


1 + Si + S2 71+1 


Si 


a _ 1 

-, S2 = - 

71 — a 71 — a 


e = S 2 /S 2 


Si TlSi 

a ’ 1 + Si 

Si _ 1 + \zi p 
aS 2 7l|z2p 


n|zip 
1 + kiP’ 
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Figure 4. P n = {(0,0), (1,0), ..., (n + 1,0), (0,1), (1,1)} 


In particular, 


and therefore 


<1 1^:1 P < 


n — 1 


TZf = <{ (zi,Z2) : 1 


n — 1 


We have 


log(l + ||e-"^/2 . zf) = log (1 + |zip + = log(l + IziH + logL±1 


n|z2p 




0,log- 


= log|z 2 p - log(l + |zip) + logn . 


Therefore 

Hence, 


\Vi + kiP’V’V"’‘''"i + l^iP. 

u = log 122 p + nlog(l + |zip) + (n + 1) log(n + 1) - nlogn . 


= n^5aiog(l + |zin 

ZTT 


^-hp{z) 


(n+ir+i iz2p(i + kii"r 


for z G TZf ■ 


n" (1 +|2:iP + k2p)"+^ 

Now suppose that z is a point in TZf'- Since Tz _L Tp', we can write Tz = (ti, uti ). Let 


q{z) = {n + • z) = (^c, i(n + 1 - c)^ G F' . 


As before, we write 

Si = kiP, S2 = |z2p. Si = |e"^i/^Zip = e'^^si, §2 = = e""^is2 . 

By (123, we have 


Si 


S2 




I + S1+S2 n+1 I + S1+S2 n \ n+1 

Solving for Si, S 2 , we obtain 


Si = 


n—1 n + 1 — c 


S 2 = 


n — 1 


(26) 


(27) 


(28) 
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Therefore, e = S 2 /S 2 = 12 ^ 2 ! ^/(n — 1), so we have 


n = -log(n- l)|z 2 | 


Si = 


l^il 


(n — 


Thus, 


|zi| 


(n-l)i/»|z2pA 


By m. 


log(l + ||e • zf) = log(l + Si + S 2 ) = log (+ 

\n — 1 


{qiz),T^) = - c)j ,(Ti,nTi)^ = (n+ l)ri = log(n - l)| 22 p • 


Hence by (1^ . 




\n—1 (n — l)i/”|z2p^’" 

(Note that ijjp has constant rank 1 on TZp/ as indicated by Theorem II. df ini. since '0p = (n + l)g*wcpi on 
TZp', where g is the multi-valued holomorphic map to CP^ given by g{zi,Z 2 ) = (c„zi, 

By Theorem ll.dr iiib we know that ipp = 0 on TZy. To complete the description of ip Pi it remains to 
describe the regions TZp' and TZy. We note that a forbidden point z lies in TZp' if and only if c > 1. By 123, 
this is equivalent to Si > or 

|z2p<(n-ir-i|zi|2-. 

Therefore 

7^F' = |(2l,Z2):^<|z2p<(n-l)"-l|zlp^ |zi|2 > ^ 

' n — 1 


This leaves us with 


7^„ = ( 21 , Z 2 ) : | 22 |^ > (n- l)"-i|zi|^”, |zi|^ > 


n — 1 


n — 1 


(29) 

(30) 


These subregions are illustrated in Figure 01 To summarize: 

(n-k 1)^95log(l -k | 2 ip -k I 22 p) = (n -k 1 )ivfs for z S Ap (see (1^ 1 


ipp = < 


n^ddlog{l + |zip) 

(n + l)^aaiOg + („_i)'l7i|,^|2/n ) 


for z G TZp 


(see (123) 
for z G TZp, (see (|23) 


0 


for z G TZZ 


(see (T^ l 


Thus, as stated in Theorem ll.3l ipp has constant rank 1 on the flow-outs of the one-dimensional faces F, F\ 
and vanishes on the flow-out of the vertex v. 


1.3. Methods of proof. Having stated our principal results, we now briefly outline some key ideas in the 
proofs. The key result is Theorem ll.2l on the mass of polynomials with Newton polytope P. To prove it, we 
begin with an easy formula 

Ei/np(I/(^)Ifs) = ^^•^p)E|AP (I/(2:)|fs) = ^j-_^p)n|Arp(2,2) , (31) 

(see @ relating expected mass to the conditional Szego kernel Hi^vp, i-e. the orthogonal projection onto 
Poly(lVP). In general, the term ‘Szego kernel’ of a space S of functions refers to the kernel for the 
orthogonal projection to S from the space of all functions; i.e., it is of the form Il{x,y) = ^ji^)sj{y)i 
where {sj} is an orthonormal basis of S. Precise asymptotics for ^^.^p^ n|jvp(z, z) are given in Theorem l4.ll 
Note that the normalizing factor is straightforward to evaluate since 

#(iVP) = dimPoly(iVP) = Vol(P)A^™ -k ... , 
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is the Ehrhart polynomial EEi. 

To analyze the conditional Szego kernels n| 7 vp( 2 , w), we introduce the polytope characters 



XNpie^'^)— Y, = (e*^U...,e*^-) . 

(32) 


aeNP 


We observe that 




n|A,p(z, w) - J Anp{z, • w)xNp{P'^) dip , 

(33) 

where 

T- = {(Cl,... Cm) G (C*)™ : 101 = 1,1 < J < m} 

(34) 


is the m-torus, and where nArp(z, w) is the Szego kernel of Poly(iVpI]), the space of all homogeneous poly¬ 
nomials of degree Np. To obtain asymptotics, we need to analyze the behavior of XNp{e'^) as iV —> oo. 

To do so, we use the Euler-MacLaurin formula of Khovanskii-Pukhlikov EE], Brion-Vergne ]WniW2] . 
and Guillemin Eu): 


Xp{e'^)=Todd{J^P,d/dh) 



{w e C™, ||w|| < e) , 


(35) 


where P{h) is of the form {a; : {uj,x) + Xj + hj > 0, 1 < j < n} (P(0) = P) and where Todd{Pp,d/dh) is 
a certain infinite order differential Todd operator. Upon dilating the polytope, one obtains 

XNp{e'") = N’^Podd{Pp,N-^d/dh) { f e^^'^’^'^dx 

\Jp(h) 

In Proposition 13.11 we show that the family xnp is a complex oscillatory integral of the form: 


< e). 


(36) 


h—0 


XArp(e“) = J e^^'^’'^'>[Ao{x,w)N^ + Ai{x,w)N^ ^ H- \-An{x,w)]dVoln{x) , (37) 

for Qw sufficiently small (but Stw arbitrary), where the Ai are analytic functions that are holomorphic in w 
and algebraic in x 

We then substitute this expression in and use the method of stationary phase for complex oscillatory 
integrals |Hol Ch. 7] to obtain the asymptotics of \d\iqp{z,z). For the case where 2 is in the classically 
allowed region, we easily find that the critical point of the phase is given by v? = 0 and x = pp.{z). Since the 
phase vanishes and has nondegenerate Hessian at the critical point, we immediately obtain an asymptotic 
expansion. The case where z is in the classically forbidden region is more subtle. Since pii{z) lies outside 
of P for this case, the phase has no critical points. To complete the analysis, we must deform the contour 
to pick up critical points. In particular, we consider the complexihcation (C*)™ of T’" and deform T™ to 
a contour of the form (log |Ci |,..., log |Cm|) = t G M™, on which the phase has a ‘critical point’ with ip = 0 
and X = Xc G dP. To be precise, the derivative tangential to the face of P containing Xc vanishes at Xc 
(while the normal derivative is nonvanishing), and furthermore the phase takes its maximal real part at Xc- 
Indeed, r is the unique vector and Xc is the unique point q{z) £ dP used in formula (1 1611 for the decay 
rate bp{z), and the maximal real part of the phase is —bp{z). We then obtain an asymptotic expansion from 
the method of stationary phase on domains with boundary. 


1.4. Brief outline and remarks. To assist the reader in navigating this article, we give a brief outline: We 
begin in Section|2by reviewing some notation and terminology involving convex polytopes and Szego kernels 
for spaces of polynomials. Section |3 states and proves our oscillatory integral formula for the characters 
Xnp{P")- The heart of the paper is Section E where we derive the diagonal asymptotics of the Szego kernel 
fTheorem 14.111 for the spaces Poly(IVP), using our character formula from SectionOjas one of the ingredients. 
(Theorem 11.21 is an immediate consequence of Theorem ED) Then applying our Szego kernel asymptotics, 
we prove Theorems imioi and o on the distribution of zeros in Section [3 We also include an appendix 
(Section EJ, where we describe an alternative geometric approach to our results using toric varieties. An 
index of notation is included at the end of the paper. 

Before embarking on the proofs, we comment on our borrowed terminology from physics and on the role 
of toric geometry: 
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1.4.1. Tunnelling of zeros. Some of the terminology we use—allowed and forbidden regions, mass density— 
is taken from the semiclassical analysis of ground states of Schrodinger operators Hn = —fi^A + V on 
R". The well-known ‘Agmon estimates’ of ground states (cf. |Ag| ) show that £^-normalized ground states 
or low-lying eigenfunctions Hfiip = Eip of H are concentrated as /i —> 0 in the classically allowed region 
Ce = {x G K" : V{x) < E} and |(/3(a;)p decays exponentially as ft ^ 0 for points x in the complement. 
In our setting, the Hamiltonian is d*d on £^-sections of powers of the hyperplane section bundle on CP”’', 
ft = 1/A^, and the ground states are the holomorphic sections. Replacing |(/3(a;)p is the expected mass density 
I/('^)IfS of the polynomials in this subspace. Theorems 11.1 iTTSI show that the Newton-polytope constraint 
on the polynomials creates a ‘tunnelling theory’ for zeros and mass. 

The decay function &p in Theorem o is analogous to the Agmon action to the allowed region |Ag| . We 
will also obtain (Proposition the integral formula 



z) +pgL{e 




• da 


(38) 


(where JJ"* denotes the integral over any path in K™ from 0 to t^), which could be interpreted as an action, 
thereby bringing the results closer to classical Agmon estimates. 


1.4.2. Role of tone varieties. Since Poly(A^P) is naturally isomorphic to the space Lp) of holomor¬ 

phic sections of a natural line bundle Lp over a toric variety Mp associated to P (see O, the reader may 
wonder whether toric varieties play any role in this paper. The answer is that neither the statements nor 
proofs of our results involve toric varieties in any essential way. However, the theory of these varieties does 
give an alternative approach to the asymptotics of the conditional Szego kernel, as will be explained 3E1 It 
was also the approach in the original version of this paper |SZ2| and motivated some of the ideas. 

The reader may also wonder how the results of this paper would change if, instead of defining the Gaussian 
measures 7 |p to be the conditional measures of one fixed ensemble, we defined the measures on Poly(A^P) ~ 
ftl°(Mp,Lp) to be the Gaussian measures 7 ^” induced by the inner product induced by a Hermitian 
metric on Lp and its curvature form comp ■ As a special case of our results in lEZH Prop. 4.4], one obtains 
(under the added assumption that Mp is smooth or maybe has orbifold singularities) 

(^/i...../ic) = ^Mp + O dij A ddj'j + O , (39) 

where Ij,dj are the action-angle variables of the moment map jjLp : Mp R.™ of the T^-action on Mp. 
Related results in a somewhat different set-up have also been obtained by Malajovich and Rojas [MR] . This 
quite different law shows that the measures 7 ^^ are singular relative to ‘^np\np in the limit as N —>■ 00 . 

As mentioned previously, the polytope P in this article is only used as a constraint on the polynomials in 
creating conditional measures. The norms of the monomials 2 ;“ are fixed (as their CP™ norms). Thus, the 
change in the distribution of zeros as P varies is due solely to the choice of which monomials occur in the 
polynomials. In the case of 7 ^^^, the norms of the monomials vary as P varies since they are £^-normalized 
on the toric variety Mp. Hence the variances of the coefficients of a polynomial in the 7 ^^” ensembles 
depend on the choice of metric on Lp. This dependence creates complicated biases towards some monomials 
and away from others as P varies, making it difficult to understand what a comparison between the 7 ^^ 
ensembles would be measuring. 

In fact, our results using conditional measures also apply to polytopes that are not convex and hence 
do not correspond to any toric variety. Indeed, suppose that P is a nonconvex ‘lattice’ polytope. Then 
Theorems II . lhn~H hold with the following modifications: The point q{z) satisfying (I14II ~(I15 (I is not always 
unique; instead we choose q{z) to minimize bp{z). Then the function bp{z) is C^, not C^, and ijjp is a positive 
current, which has singular support. Part (iii) of Theorem ll. 31 applies only to the absolutely continuous part 
oi ^pp. More significantly, the limit measure in Theorem ll.ll does not vanish on the forbidden region; instead, 
it is orthogonal to volume measure there. This observation is relevant to the theory of ‘fewnomial’ systems, 
which we discuss in a forthcoming paper |SZ4| . 


1.4.3. Acknowledgements. Our interest in polynomials with a fixed Newton polytope was in part stimulated 
by a discussion with A. Varchenko at the outset of this work. We would like to thank M. Brion for many 
helpful comments regarding polytopes and the Euler-MacLaurin formula. We also grateful to T. Theobald 
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for giving us permission to use Figure 0 from his paper and A. Carass for providing the illustrations in 

Figures 000 El and El 


2. Background 

2.1. Fans. By a convex integral polytope, we mean the convex hull in ffi.™ of a finite set in the lattice Z™. 
A convex integral polytope P with nonempty interior P° can be defined by linear equations 

:= {x,Uj) + \j>0, {j = l,...,d) , (40) 

where Uj € Z"* is the primitive inward-pointing normal to the j-th facet (codimension-one face) 

:= {x G P : ijix) = 0, £k{x) > 0 for k ^ j} . 

For each point x G P, we consider the normal cone to P at x, 

Cx = Cf := {m G K™ : {u, x) = sup(M, y)} , (41) 

y&P 

which is a closed convex polyhedral cone. We decompose P into a finite union of faces, each face being an 
equivalence class under the equivalence relation x ^ y Cx = Cy. For each face F, we let Cp denote 

the normal cone of the points of F. Note that by our convention, the faces are disjoint sets. We 
shall use the term closed face to refer to the closure of a face of P. 

Each face of dimension r (0 < r < m) is an open polytope in an r-plane in R™; i.e., the 0-dimensional 
faces are the vertices of P, the 1-dimensional faces are the edges with their end points removed, and so forth. 
The facets Fj = Fp~^ and their normal cones are given by: 

Fj = {x G P : £j{x) = 0} , Cpj = {—tuj : f > 0} . 

The TO-dimensional face is the interior P° of the polytope with normal cone Cpo = {0}. 

For each x G P, we let 

f7(x) = {j G Z : 1 < j < d, £j{x) = 0} . (42) 

One easily sees that Cx = Cy J{x) = J {y), and hence we can write J (P) = J{x), where x lies in 

the face F. In particular, J{FJ^~^) = {j} and J{P°) = 0 . The polytope P is called simple if ffj{v) = m 
for each vertex v of P. In this case, {uj : j G Jiv)} is a basis for R.™ for each vertex v (since P° 7^ 0 ), and 
furthermore ffJ{F) = codim P for all faces P. The polytope P is said to be Delzant if {uj : j G J{v)} 
generates the lattice Z™ for each vertex v of P. 

The convex integral polytope P determines the fan Fp := {Cp : P is a face of P}. A fan F in R™ is a 
collection of closed convex rational polyhedral cones such that a closed face of a cone in F is an element of 
F and the intersection of two cones in P is a closed face of each of them. (Fans are used in the algebraic 
construction of toric varieties; see [Fill §1.4].) An example of a convex integral polytope and its fan is given 
in Figure 0 



Figure 5. A convex polytope and its fan 
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We shall also consider convex integral polytopes P with empty interior. In this case, the faces of P and 
the normal cones are defined exactly as above. If P° = 0, then the normal cones of P all contain the 
linear subspace of R"* orthogonal to P (contrary to the case P° 0, where the normal cones are pointed, 
i.e., do not contain any lines in R"*). A convex polytope in R™ of dimension n < to is also said to be simple 
if it is linearly isomorphic to a simple polytope in R" or equivalently, if every vertex is the intersection of 
precisely n closed faces of dimension n — 1. 

2.2. Szego kernels. By homogenizing the polynomials in Poly(P), we obtain a hnite dimensional subspace 
S of £^(5’^’"+^). By the ‘Szego projector,’ we mean the orthogonal projection : £^(A') ^ S. We recall 
that the projector II^ is given by a kernel of the form 

kN _ 

n5(a::, v) (^)sf (y) ’ (43) 

i=i 

where {sy} is an orthonormal basis of S. 

We now describe two sequences of Szego kernels, which play a crucial role in our main results: 

2.2.1. The projective Szego kernels. As a first example, we recall that Poly(pS) can be identihed with the 
space of degree-p homogeneous polynomials in to + 1 variables by identifying the (non-homogeneous) poly¬ 
nomial 


f {Zl, ■ ■ ■ , Zm) — ^ ( 

|a|<p 


CaZ 


iZ = Z, 


■■Z„ 


0 


with the homogeneous polynomial 


F{Co,-..Xm)= CaCo”'“'Cr---C 

l«l<p 

We equip the space Poly(pE) with the inner product on 


{f,9) = 


FGdiz = 


ml Jg2m+1 


m\ Jc 


where 


wfs = 


f{z) 9 {z) 

c- (1 + lkP)^ 

27T 


‘^Fs(P. />5 e Poly(pE), 


5aiog(l -f ||zp) 


(44) 


(45) 


is the Fubini-Study Kahler form on C™ C CP™. 

A basis for Poly(pS) consists of the monomials Xa(z) = • • ■ z"^, |q!| < p. The monomials {xa} are 

orthogonal but not normalized. Their £^ norms given by the inner product 14411 are: 


IIXcll = 


(p - |a|)!ai! • • • Om! 


(46) 


(p -I- to)! 

Note that the norm || • || depends on p. Thus we have an orthonormal basis for Poly(pE) given by the 
monomials 


1 


IIXo 


■Xa = 


(p + to)! 


where 


(p - |a|)!Q;i! • • -ami 

'P 


{p + my. (p'' 

Xa = \ -i- Xa , 


p! \a 


|a| < p . 


p! 


\C^J (p - |a|)!Q:i! • • • Om! 

We let : ^^m+i ^ denote the homogenization of Xa'- 

Xa(a:) = •• - 0;“- . 

Hence the Szego kernel Hp for the orthogonal projection to Poly(pE) is given by: 

^p{x,y) = ^ ||^^||2 Xa(a:)xa(p) = {x,yY , 


(47) 

(48) 

(49) 


|a|<p 


IIXo 


(50) 
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for x,y G (The sum is the usual Szego kernel for the sphere; see |BSZI §1.3.1].) 

We also identify the point z G (C*)™ with the lift x = (1; zi,..., Zm) € 5 ' 2 m+i^ write 


XS(^) = 

The Szego kernel can then be written explicitly as 


(1 


np(z, w) = 


{p + m)\ 


2)p/2 • 


{p + m)\ 


p\ 


(1 


2)P/2(X 


||w||2)p/2 


1 + {z,w) 


p\ [(l+|kP)l/2(l+||t^||2)l/2_ 


(51) 


(52) 


2.2.2. The conditional Szego kernels associated to a polytope P. In this case, the relevant space of polynomials 
is the subspace Poly(P) C Poly(pE) of polynomials with Newton polytope P. Here, we may choose any 
p > degP := max{|a| : a G P}, but we normally choose p = degP. In the conditional Szego kernel, we 
retain the Fubini-Study inner product on this subspace. Hence this example is very similar to the previous 
one. The main difference is that an orthonormal basis of Poly(P) is given by 

Definition: The conditional Szego kernel Hjp is the kernel for the orthogonal projection to Poly(P) with 
respect to the induced Fubini-Study inner product: 

1 


H 


\P 


ix,y) = 


aeP 


IXo 


TXaix)xUy) ■ 


(53) 


When defining the term ‘random polynomial with fixed Newton polytope P\ we wish to use an >C^-norm 
on monomials which is defined independently of P. This explains why the conditional Szego kernel is the 
essential one in our problem. The conditional Szego kernel can be written explicitly on C"* as 

iP + my. T,aGP' 


n|p(z,w) = 




p! (l + ||z| 


2)P/2(X . 


W 


2)p/2 


(54) 


3. Polytope character 

This section is devoted to the asymptotic analysis of the polytope characters (P|l. Our interest is in the 
asymptotics of the ‘ray’ of characters Xivp(e“) defined on (C*)™ by 

y^p(e“) = Y w G C™. (55) 

aGNP 

We shall derive a formula expressing the character XNp{d^) as an oscillatory integral over the original 
polytope P with the same phase (Wjx) as in formula (ESj) of iK’FiiwninjTi . which is the starting point of 
our derivation. Although (IS31l holds only for small w, our formula holds by analytic extension for all re in a 
‘strip’ of the form 

S(e) := {it G C™ : \^Wj \ < e for 1 < j < m} . (56) 

Our integral formula for xnp{oP) is given by the following proposition, which is also of independent 
interest. 

Proposition 3.1. Let P be a simple integral polytope in R™ of dimension n (1 < n < m). Then there 
exists e > 0 such that the characters xnp can be given as integrals over the polytope P of the form 

XNp{e^) = j + Ai{x,w)N'^~^ + ■ ■ ■ + An(x,w)]dYo\n(x) , for all w € S{e) , 

where the Ai are analytic functions on P x S(e) that are holomorphic in w and algebraic in x, and 

i) Ai{x,w) G K whenever w G R.™, 

ii) Ao(a:,0) = 1, 

iii) Ao(t,u>) G whenever w is in the normal eone to P at x, i.e., whenever w G Cx C M™. 
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As discussed in the introduction, Pronosition id.l I wili later be used to obtain asymptotics of the conditional 
Szego kernel through formula When w = 0, the character fKK|i equals the number of lattice points of 

NP and Pronosition Id. 1 1 reduces to the Ehrhart formula m-- 

n 

#{NP) = ajN^-^ , n = dim P , uq = Vol„(P) . (57) 

3=0 

We note that one can easily extend to an arbitrary convex integral polytope P (not necessarily simple) 
by triangulating P into n-simplices {Ai,..., A^} whose vertices are vertices of P (see [GKZI pp. 214-216] 
for an elementary proof of this fact), and then applying the inclusion-exclusion principle. 

Remark: The characters Xivp(e“) have a natural interpretation in toric geometry, and xap(1) = ^{NP) 
can also be given by the Riemann-Roch formula on the toric variety associated to the polytope P. See 
for further details and references. 


3.1. Proof of Proposition l37ll We first reduce to the case where P has nonempty interior: If dimP = 
n < m, then by a translation we may assume that 0 S P. Next by an orthogonal projection to the n-plane 
(P) spanned by P, we may assume that w G (P). Since (P) n Z™ « Z", we may assume that (P) = R". 
Hence we may assume that n = m. 

Recalling (Uni) . we express P in the form 

P = {x : {x, Uj) -f Aj > 0, j = 1,..., d} , (58) 

where the vectors Uj are the inward-pointing normal vectors to the facets of P. The normals Uj are normalized 
so that Uj is a primitive element of Z™ (i.e., ruj G Z™ r G 1^). 

We shall use the formula given independently by Brion and Vergne |HV1 1 Theorem 3.12] and by Guillemin 
|Gu[ Theorem 4.1] : 


Xp{e^)=Todd{P,d/dh) 



h—0 


wGC^ , Ur;]] < ep. 


(59) 


Here, 

P{h) = {x : {uj,x) + Xj + hj > 0, 1 < j < d} , (60) 

P = Pp is the fan associated to P, and ep is a positive constant. (E.g., if P is Delzant, then ep = 27r.) Also, 
Todd(iF, d/dh) is the (generalized) Todd operator constructed by Brion and Vergne |Wn pp. 374-376] as 
follows: 

Eor each cone Cp G P, we consider the open fc-parallelogram 


Up : 0 < < l| , 

where k = codimP = dimC'p’. (Recall (03).) We let U denote the (disjoint) union of the Up, where P 
runs over the faces of P (including the open face P°, where Upo = {0}). Equivalently, P is a ‘bouquet of 
m-par allelograms ’: 

U = [J : 0 < < l| , 

vertices v 

and is an open neighborhood of 0. We consider the finite set 


:= p n z™ . 


(It is easy to see that Tp = {0} if and only if the polytope P is Delzant.) For each lattice point j GTp, 
we define a unique vector 5 ( 7 ) = ( 31 ( 7 ),... , 5 m( 7 )) G ([0,1) n Q) as follows: Let P be the face such that 
7 G Up; then 

m 

l = Yd3il)'^U 5j(7)7^0 
1=1 


J G J{F) ■ 
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We note that the map 7 1 -^ 5 ( 7 ) is injective. The Todd operator is then given by: 

d 

Todd(J^, d/dh) = E n 

7Gr^ 7=1 

As noted in mu (Remark 2.16), if P is a Delzant polytope, then Todd(iF, i9/9/i) is the usual Todd 
operator 

Todd(9/a/i) = W^^^Toddid/dhj), Todd(z) = Todd(l,z) = z(l - (62) 

and formula H59|l is due to Khovanskii and Pukhlikov |KP] for that case. (Formula II59II has a generaliza¬ 
tion to arbitrary convex rational polytopes; see m^-) Since the term in m corresponding to 7 = 0 is 
the usual Todd operator and the other terms vanish at the origin, the constant term in the series expan¬ 
sion of Todd(iF, d/dh) equals 1 for all simple polytopes. Furthermore, the coefficients of the expansion of 
Todd(P, d/dh) are real since if we regard the g{'j) as elements of then the set { 5 ( 7 )} is a union of 

subgroups and hence is invariant under multiplication by — 1 . 

Let us check how formula dilates. The fan does not change when the polytope is dilated, so the Todd 
operator is the same for P and for NP. Noting that 

{NP){Nh) = {x : {uj,x) P NXj + Nhj > 0} = N{P{h)) , (63) 

we change variables h 1-^ Nh, x Nx in (El to obtain 

XNp{e^) = Todd{P,N-^d/dh) ( f 

\Jpih) 


|w|| < ejr 


(64) 


h—0 


Todd(a,z) = ^ . (61) 


To prove Prnnosition 13.11 we use the following lemma to change variables so that the integral is over the 
fixed polytope P while the integrand depends on h. 

We say that a map f : P ^ R"* is a regular rational map if it can be written in the form / = 
(Pi/^ii ■ ■ • ^Pm/gm), where the pj and qj are polynomials on R"* and the qj are nonvanishing on P. 

Lemma 3.2. Let P be a simple polytope in R*” with nonempty interior. Then there exist regular rational 
maps fj:P—> R™, 1 < j < d, such that for h G R'^ suffieiently small, the map jh ■ P R™ given by 

d 

lh{x) = x + J2hjfj{x) 

1=1 

is a diffeomorphism from P onto P{h) mapping each face of P onto the corresponding face of P{h). 


Proof. Recall that the assumption that P is simple means that there are precisely m hyperplanes {x : 
(x, Uj) + Xj = 0} through each vertex, and hence for small h, the faces of P{h) correspond to the faces of P. 

We first construct an ‘algebraic partition of unity’ {ipk} over P as follows. Let ..., n'*} denote the set 
of vertices of P. Recalling innii . we let tpk ■ P ^ [0, -boo) be given by 

=Y[{^j ■■ j ^ J{v'')} , k = l,...,s. 


We then let 


i^k = 


ipk 


We easily see that 0 < ■i/'fe < 1, Ek=i i^k = ^ and that i/’^T^(0) is the union of the closed facets of P that do 
not contain v^. 

The vertices v^{h) of P{h) are of the form 

v'^{h) = + L{k)h , 1 < fc < s, ||d|| < £ , 


where L{k) : R”* ^ R™ is a linear transformation. (In fact, it is easy to see that L{k) is given by first 
projecting onto the m coordinates corresponding to the facets incident to Vk, and then applying an element 
of GL(to,R).) Let 


= L{k)E^ G R™, 
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where {E^} is the standard basis of so that 

d 

v^{h) =v^-\- ^ hjW^^ , ||/i|| < e . (65) 

For each k, let be the vertices connected to by an edge (1-dimensional face) of P. 

For 1 < fc < s, 1 < j < d, let fj-j : ^ R™ be the unique affine map such that 

fkjiy'-) = for l = k, v{k,l), v{k,m) . 

So by lIHKll . if Z = fc or if is a vertex connected to by an edge of P, then 

d 

v'-{h) = v'- + '^hjfkj{v’‘) , for I = k, i^(/c,l), v{k,m) , ||/i|| < e . ( 66 ) 

i=i 

We claim that the maps 

S 

fj ■= ' 0 fc/fei 
fc=i 

satisfy the conclusions of the lemma. To verify the claim, we let 

d 

Ihix) =x + ^hjfkjix) , l<fc<s, 

so that 

lh = Y^ il’kll ■ 

It follows from linni) that when a: is a vertex we have 
since ipi{v^) = 0 when k ^ 1 . 

Next we consider points on the 1-skeleton of P. Suppose that x = tv^ -I- (1 — t)v^, where 0 < t < 1 
and are joined by an edge of P. Then tphix) = 0 for k > 2 and we again conclude from 16611 that 

7ft and 7 ^(u^) = v‘^{h) for k = 1,2. Since the 7 ^ are affine maps, it follows that 

7 ft(a;) = tv^{h) -I- (1 - t)v‘^{h) 

so that 7 ft takes each edge of P linearly onto the corresponding edge of P(h). 

In general we do not have linearity on faces of dimension > 1, so we need to modify the argument for 
the higher skeletons. Suppose now that x lies in a 2-dimensional face F with vertices v^,... ,v'' (t > 3) and 
X = tiv^ + t 2 V^ + t^v^, where mintj > 0 and ^ tj = 1. Let 1 < / < r. Since the points v^{h ),..., v'^{h) lie 
on a plane St and the 7 ^ are linear, it follows from the above argument that 7 ^( 1 ^*) S Sh for 1 < fc < r. (For 
example, if are each connected to by an edge, then 7 ^(u*) = v\h) for I = 1,2, 3, but for 3 < / < r, 
we can conclude only that 7^(w0 ^i^s on the plane through v^{h), v‘^{h), v^{h).) Since for k > r, there exists 
a closed facet containing vi,V 2 , V 3 but not containing Vk, and hence ipk{x) = 0. Therefore, 

7ft(a^) = ^3'^k{x)-fl{v^) & Sh ■ 

1< j,A;<r 

Furthermore, for h small, D^h ~ D^q = I, so "fh\F '■ F Sh is a. local homeomorphism. Therefore, 
9Shilh{F)) = jhidF), which is the boundary of the convex polygon F{h) with vertices v^{h), ...,v'"{h). It 
follows that ^h(,F) = F(h) and hence 7 ft|-p : F F{h) is a (global) homeomorphism for sufficiently small h. 

Continuing by induction on dim F, we conclude that ^h '■ F k, F(h) for all faces F of P, and in particular, 
7 ft maps P diffeomorphically onto P{h), for h small. □ 

Completion of the proof of Provosition l,V. 1 \ Let 'jh ■ P ^ P{h) be as in Lemma |3.2I Making the change of 
variables x 1 -^ lh(x), we have 

f e^^^’^Ux= f detDjhix)dx = f h) dx , 

Jp(h) Jp Jp 


( 67 ) 
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where G is a regular rational function on P x K of the form 

G{x,h)= Go,{x)h 

Hence by (|P| and we have 

XNpie'^) = J Todd(P, lV-i5/5/i) (e 


Gn = 1 . 


N[{w,x)+^j{w,fj{x))hj] 


G{x,h)) 


h=0 


dx , 


\w\\ < e . 


( 68 ) 


To justify interchanging integration and Todd differentiation in the above and to simplify (PI , we use 
the following identity from [KPI Lemma 1]): Suppose that G{h) G , hk] is a polynomial and F is a 

convergent power series about 0 G with domain of convergence Then 


F{d/dh) =G{d/dq) 


qefl, heC^. 


(69) 


(Equation llH!7ll is easily verified by noting that both sides equal F {d/dh) G (d/dq) exp(S^^^( 7 j 7 ij).) 
We apply lIMIl with F{q) = Todd(lF, A^“^g), G{q) = G{x,q) to the integrand of ll^ : 


Todd{F,N-^d/dh) 


^N[{w,x)+i:j {wjj {x))h. 


^G{x,h) 


h=0 


= Todd{F,N-^d/dh) fG(a;, 

= e^^'"’^^G{x,d/dq) [Todd(F,iV-ig)e^i=i«^'‘^ 

G{x,N-^d/dq) Todd(F,g)|. 


h^O, qj^N{wJj(x)) 
h^O, qj^N{w,fj{x)) 


= e 

= e 


(70) 


In particular, if we let Todd^ denote the sum of the terms of the Todd series of degree < k, then 
Toddfe(F,iV-i5/9/i) [e^[{-^<^)+^j{-^,hG))hdG{x,h)] 


= e^^^’^'>G{x,N-'^d/dq)ToddkiJ^,q)\ 


Qj = {^JjG)) 


ToddiF,N-^d/dh) [em{^,^)+^A^Jjio^))hdG{x,h)]\^^^ 

uniformly for x G P and w sufficiently small, which justifies the interchange of integration and Todd differ¬ 
entiation in PI - 

Let A(a;) denote the d x m matrix whose rows are the fj; i.e., A(a;) : ^ is the linear map given by 

{A{x)y)j = {yjj{x)). By IZOI 


Todd(F, N-^d/dh) L^[i^,-)+^F^JiG))hdG{x, h)] = V N-‘Ti{x, A{x)w) , (71) 

J h=0 


1=0 


where Ti{x,q) is a convergent power series in q with coefficients that are regular rational functions on P. 
More specifically. 


To(a;, 9 ) = Todd(F,g) , Ti{x,q) = Gi{x,d/dq)Todd{F,q) {I < I < m) , 


where Gi{x,d/dq) is a homogeneous differential operator of order I whose coefficients are regular rational 
functions of a; G P. 

It follows from (EHI and II71II that the integral formula of Proposition Id. II holds for small w, with 


Ai{x,w) = Ti{x,A{x)w) . (72) 

Moreover, from the construction of the P;, we see that the integrand has no poles if ||A(x)3w|| < ep for all 
X G P. Thus by analytic continuation, the identity holds on S{e), where e = sp/ sup^^p ||A(a:)||. 

Statements (i) and (ii) follow from IZ3) and the fact that the coefficients of the Todd function are real 
and the constant term is 1. 

To verify (iii), we fix a; G P and w G Cx, and we let yj = (w, fj{x)), for j = 1,... ,d. 

Claim 1: 
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i) Vj > 0, for j = 1,..., d; 

ii) yj = 0, for j ^ J{x). 
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Proof of Claim 1: (i) Let F be the face containing x. Then x — tfj{x) is on the face F{h) corresponding to 
F of the polytope P{h), where hj = —t (for f > 0 sufficiently small) and = 0 for fc ^ j. In particular, 
x — tfj{x) G P. Since w is in the normal cone to F, we have {w,x — tfj(x)) < {w, x)\ i.e., yj = {w, fj{x)) > 0. 

(ii) If j ^ J{x), then F{h) C F for h as above, and hence fj{x) is tangent to F, so that yj = {w, fj{x)) = 
0. □ 

Claim 2: Let Gp = D Up, where F is the face containing x. Then 

Aoix,w)=J2 n Todd(F2"9^(^),y,). 

76Gf iGj'(F) 


Proof of Claim 2: We have 

Ao(a:,w) = Todd(F,?/) , y = (yi,..., j/^) , yj = {w, fj{x)) . 

We first note that 

Gp = {-fGZ'^:-f= tjUj, 0 <tj < 1 } . 

3&J{F) 

Thus, if 7 e TjP \ Gf, then there exists j ^ J{F) such that 5 ^( 7 ) yf 0. But by Claim 1, the corresponding 
yj = 0. Hence the term corresponding to 7 in the sum (I till! for Todd(F, y) vanishes. 

Now suppose that 7 G T^ C Gp. Then by definition, 5 ^( 7 ) = 0, for all j ^ J{F). Therefore 

= Todd(l, 0) = 1 for j ^ J{F) . 

Hence the sum (HU for Todd(F, y) reduces to the sum in Claim 2. □ 

We easily see that the set 

Qf ■■= {{9j{l))jeJiF) e Q7Z" : 7 e Gf} 

is a subgroup of Q”/Z"', where n = ffJ{F) = codim F. The final conclusion (hi) of the Proposition is then 
a consequence of Claims 1-2 and the following elementary lemma. 

Lemma 3.3. Let L 6 e a finite subgroup of T” and let yj > 0, for 1 < j < n. Then 

n 

Y IlTodd(aj,yj) > 0 , a = (oi,..., a„) . 

oGF j=l 

Proof. We can assume without loss of generality that the yj are positive. For suppose that y„ = 0, for 
example. Then Todd(a„,y„) = 1 if a„ = 1, and Todd(o„,y„) = 0 if a„ yf 1. Then the sum reduces to the 
corresponding sum over the group V ■.= {a’ G : (o', 1) € T}. 

Let rj = e~y^ < 1. Then 

n n 

^ ■= Y. n Todd(a„y,) = ^ f[il - a,r,)-i = ^ ^ = Y > 

aGF j=l qGF 3 = 1 aGF /3 gN" /3GN" 

where 

cp = Y^^ ■ 

aGF 

We claim that cp is either 0 or o(r). Indeed, if = I for all a G F, then cp = o(r). On the other hand, 
if there exists 6 G F with yf 1, then cp = — SaGF {b ■ a)^ = bl^cp, and hence cy = 0. Since 

Co = o(r), it then follows that S' > 0. □ 
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4. Mass asymptotics: proof of Theorem II .21 


In this section, we prove a precise asymptotic formula for the conditional Szego kernel on the diagonal 
1 Theorem 14. Ill , which yields the mass asymptotics of Theorem 1 1.21 

First, we verify formula m for the expected mass density. Let P denote a convex integral polytope 
in R™. Recalling the definition 0 of the conditional probability measure 7 p|p on the space Poly(P) of 
polynomials with Newton polytope P, we see that the expected value of the mass density with respect to 
7 p|p is given by: 


E|p 


(i/(^)Ifs) — 

a,l3GP 


^iK>^l3)Xai.z)xi3{z) 

IIXo||k/3||(l+||zP)P/2 ■ 


Since the Aq, are independent complex random variables with variance 1 (i.e., E|p(AqA^) = df), we have by 
(I5d|l : 

E|P (|/(^)|ls) = E • (73) 

It then follows by expressing the Gaussian in spherical coordinates that 


E; 


ep(I/(^)Ifs) — (I/(^)Ifs) — Trp^T(^>^) ■ 


1 


Replacing P with NP, we obtain 

E.„p(|/Wlls) = ^^n|^p(z,z). (74) 

Recall that by the number ^{NP) of lattice points in the polytope NP is a polynomial in N of degree 
equal to the dimension of P. The mass asymptotics of Theorem o is an immediate consequence of (El) 
and the asymptotic expansion of the conditional Szego kernel on the diagonal given in Theorem 14.II below. 


Theorem 4.1. Suppose that P is a convex integral polytope in such that P C pE and P (f_ 5(pE). Then: 

i) If P° ^ 0, then for z in the classically allowed region Ap, we have 

m 

IllNp{z,z) = Y[{Np +j) + Rn{z) , \\RN\\c>‘iK) = /or fc = 1, 2, 3,..., 

i=i 

for all compact K C Alp, where Xk > 0. 

ii) On each open forbidden region TZp, 

n|^p(z, z) = K(^) + c({z)N-^ + • • • + cf {z)N-’^ + Pf (z)], 

where r = dim F and 

(a) c^ G C°°{TZp) and Cq > 0 on Ti°p; 

(b) \\Rk\\ci(K) = for all compact K C TZp and for all j, k; 

(c) bp>0 on{C*)'^\A^; 

(d) bp is given by formula \lt)\) : 

(e) bp G Cg((C*)"*) (with bp = 0 on Ap), and bp is C°° on each closed region TZp. 

Note that bp fails to be at transition points, since the limit expected zero current ifp of Theorem ll.dl is 
discontinuous at the transition points and is given in terms of the second derivatives of bp by formula ESI- 


Remark: For the case where P = {/3} is a single lattice point in pE°, the allowed region is all of 
(C*)™. Since the point q{z) given by (IHll lies in dP, in this case g(z) = /3 for all z G (C*)"*. Recalling 
we then have 


Theorem o then says that 
^\N{0}{z,z) = 


=-21og7VJ/3(z) = log x^(/r -logx))(z) 


IIxpK 


:\X0iz)\ps = N ^ {co + ciN +C 2 N H- )M0(z) 


\2N 


(75) 

(76) 
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(The asymptotic formula 117611 is also an immediate consequence of formula 14611 for the C? norms of monomials 
and Stirling’s formula. A more precise formula and general results on the decay of monomials normalized 
over toric varieties are given in |STZ2j .') 

Replacing /3 with an arbitrary (non-lattice) point a; G P° in 17511 . we define 

:=-21og7W,(z) . (77) 

Recalling m, we see that for a convex lattice polytope P, we have 

6p(z) = 6g(^)(z) , n|jvp(2,2) = -hcf(0)A^"^ H-] . 

To compute fe{a;}, we use (P|l with q{z) = X and /i(e • z) = ^x, which yields 

Tj = log \zj\^ - logXj + log(p - Xj) , 

and hence 

AIL r \z\^^ 

log - log ^ ||^||2)p (^o=p-Er=i^j) • (78) 

The asymptotics of Theorem EH are uniform away from the transition points only. To take care of the 
transition points, we shall also prove the following local uniform convergence result on all of (C*)"*: 

Proposition 4.2. Let P he as in Theorem \4.1\ Then 

^ logn|wp(z, z) -hp{z) 

uniformly on all compact subsets o/(C*)™. 

Proposition EH is the main ingredient in our proofs of Theorems [QHn on the asymptotics of zeros in 

dHl 

The remainder of this section is devoted to the proofs of Theorem 14.1 l and Proposition^^! To motivate the 
argument, we begin by using our integral formula for the polytope character (Proposition |OJ to show that 
the integral is rapidly decaying when z is in the classically forbidden region, while for z in the allowed region, 
the phase is of positive type with a nondegenerate critical point and hence Ii\p!p{z, z) has an asymptotic 
expansion in powers of To obtain our expansion in the forbidden region, we use Pronosition llTlI to give an 
oscillatory integral formula for the conditional Szego kernel n| 7 vp(^:, z). We then seek deformations 

of the contour of integration so that the (analytic continuation of) the phase picks up a critical point. In 
an we formulate necessary and sufficient geometric conditions for the existence of a critical point where 
the phase has maximal real part, and in ao:Ti we show that there is a contour for which these conditions 
are satisfied. Assuming that P is simple, we obtain in H4.2.2I the asymptotic expansion (ii) by performing 
a stationary phase analysis over the polytope when the critical points of the phase lie on the boundary. 
In we show that the decay function satishes properties (c) and (e). In rmi we verify the precise 

asymptotic formula (i), and in we complete the proof of part (ii) of Theorem EH by extending our 

results to non-simple polytopes. We prove Proposition 14. 21 in 

We start our proof of Theorem EH with a simple integral formula for the conditional Szego kernel: 

n|Afp(^:,^:) = J IlNp{z,e^‘^ ■ z)xNp{e''‘^)dip , (79) 

where we recall that 

T™ = {(Cl,... Cm) e (C*)"* : 101 = 1, 1 < J < m}. 

The identity lO follows immediately from the explicit formulas (1^ and (1^ for the Szego kernels and the 
definition itCT) of the character xnp- 

Substituting formula (1^ for Ii]s[p{z, • z) in (TTHll . we obtain: 

{Np-\-m)\ f 


-I I \-^7Tl 

1 + E,,=i e 


-iipj 


Zi 


1 + lkll 


n|7vp(^,^:) 


XNp{e"^) dp . 


( 80 ) 
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We now assume that P is a simple polytope. In H4.2.5I we shall derive the general case as a corollary 
of our results for simple polytopes. 

4.1. The classically allowed region. In order to motivate our argument, we give a short proof of part (i) 
of Theorem lO with the exponentially small error term replaced with an 0{N °°) error. The actual proof 
of part (i) will be given in H4.2.4I and does not depend on the results of this section. 

Assume for simplicity that P° ^ 0. We first choose a cut-off function ^ G C°°(R”^) vanishing outside the 
e:F-ball and with ^ = 1 on the (e;F/2)-ball. Since nArp(z, • z) decays exponentially (in TV) when ||(/j|| > Ejf, 
it follows from Pronosition 1,4. II and (I8bll that 

n|jvp(z, z) ^ J J H- ]dxd(p, (81) 

where a{N) = (p + ■ ■ ■ (p + ^) = + aiN~^ -I- • • • -I- Aq{x, 0) = 1 and the phase dbt is given 

by 


'^a{.V,x]z) = i{(p,x) -bplog 


'l + Er=ie-“' 


1 + lkll 


(82) 


Here, means modulo 0{N °°). 

The integral in m is an ‘oscillatory integral’ with ‘complex phase’ to evaluate it by the method of 
stationary phase IHol Ch. 7], we must find the critical points of dbi. The (interior) critical point equations 


dx'^A = = 0 , d;^'l’A\ip=o = ix- 


ip 


1 + lkll 


r(|zip,..., \zm\‘^) =ix- ipn^iz) = 0 


yield 

(p = 0, X = pfj.{z) . 

We note that 3 (?'I '_4 < 0 and 'I '_4 = 0 at the critical point. 

The Hessian Tidbi of (with respect to the variables p, x) is of the form 

/ C il 

n'i>A\io,x) = 

\il 0 

where I denotes the m x m identity matrix. The matrix C is given by 


(83) 


(84) 


Cjk — 


d^^A 


dipjdifk 


{0,x) 


= p{-I,6^ + I,Ik) 


where 


1 + IUIP 


; i.e., /^(z) — (Ji,..., Im) ■ 


Hence 

and the inverse Hessian is given by 


det7f«'A|(0.r>;) = 1 


(85) 





Hence, the Hessian operator is given by 

H = —i ^ d^/d(pjdxj + ^ Cjkd^/dxjdxk- 

3 jk 

By (IHU and |Hbl Theorem 7.7.5], we have a stationary phase expansion of the form 

1 


n|^p(z,z)^p™a(iV) ^ 


2k>3j>0 


j\k\2^ 


pfra+o-k jpkgJ^ 


ip—0,x—p^{z) 
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where gs is the third order Taylor remainder of the phase We note that is a function of cp alone, but 
that every tp derivative is accompanied by an x-derivative in H. It follows that only the very first term of 
the summation is non-zero. Thus, the stationary phase expansion is simply p'^a{N)N^. 

4.2. The classically forbidden region. By (CTl . dt has a critical point if and only if z is in the classically 
allowed region. Thus, in the classically forbidden region, there are no critical points and the integral m is 
rapidly decaying. 

We now let z be a point in the classically forbidden region. To evaluate the integral by the method of 
stationary phase, we need to deform the contour of integration to pick up a critical point with maximal real 
part along the contour. We complexify the real torus T™ to (C*)"* with variables Q = so that 

may be written as 

{Np + m)\ f + SjLi Cj ^ /an TT 


^Np{z,z) = 


{Np)l{2TTi) 


i+Er^iCri^.i 

i + iNip 


xNP{c) n 7 ^ ■ 

,=i 


( 86 ) 


Since the integrand is holomorphic in ^ G (C*)"*, we can deform the contours in and instead integrate 
over a torus of the form 

T™ := {C G iC*r : ICil = - , IC-I = e"*"} • 

We shall show below how to choose the parameter r = (ri,... ,rm) G ffi.™ to obtain our required critical 
point. 

Deforming lIHHIl to T™, we have 

{Np + m )! 


^Npiz,z) = 


{Np)\{2TTy' 




where 


■^o{tp;T,z) =p\og 


l + ||z||2 


(87) 


(88) 


Using a partition of unity {^ 1 ,^ 2 } of T"*, where ^i(</5) has support in a small neighborhood of 0 and is 
identically 1 near 0, we decompose 1H3 into two integrals: 


TT / N {Np + m)y^^ 

n|Arp(^, zj — j 


•=_ 

" ■ (27r)- 




(89) 


We shall show that has the desired asymptotic expansion (when z is not a transition point) and that 
= 0(e“'^X^) for some ^ > 0. 

Let dimP = n. Under our assumption that P is simple. Proposition 10 lets us write the primary part 
Xy of the Szego kernel given in (18911 as an oscillatory integral over T™ x P: 

iy = N^J j e^'^^^’^'^'^'>A{p,x,N-,T)dxdp , 

where the phase is given by 

f 1 ■ 

^{p,X]T,z) = {t + ip,x) +p\og \ — 

and the amplitude is given by 

A{p,x,N-,t) = ^i{p)[Ao{x,T + ip) + Ai{x,t + i 
We therefore look for complex critical points of 4'. We note that as before, by the triangle inequality, 

3?'I'((^, x; T, z) < 5i'I'(0, x; T, z) , for p^O. (93) 

Thus for each fixed r and z, the maximal value of Jidt occurs where p = 0. From iPmi . we obtain 


c, N;t) dx dp , 

(90) 



0-1 1-5 1 

1 + lkP ) ’ 

(91) 

-1- • • • -b An{x, t + ip)N~'^] . 

(92) 


IT- *7* 

a„W = IX — - -=-^ 

l-(-^ e ’■3 P 




-Ti-Hpi 


Zil 


— Tm — 


■ 


(94) 










RANDOM POLYNOMIALS WITH PRESCRIBED NEWTON POLYTOPE 


27 


Hence the critical-point equation a;; r, z) = 0 is equivalent to: 

- X = ix(e-^■ z) . (95) 

V 

From the second critical point equation 

dx'^ = T + iip = 0 , (96) 

it follows that there are no critical points on the interior of P, since by we must choose r 7^ 0 to obtain 
a critical point when z is in the classically forbidden region. 

Suppose r S K"* and G P. Recall that r is in the normal cone to P at x^ if (r, a;°) > (r, y) for all y G P. 
Since 'I'(t, 0, a;, z) = (r, x) + Cr,z, where Cr^z is independent of x, it follows that takes its maximum along 
the contour at the point {v? = 0, x = a:°} if (and only if) r is in the normal cone C^q to P at a;°. Recall 
that if T is in the normal cone at x^, then r must be orthogonal to the face of P containing x^. 

Therefore, to obtain a ‘critical point’ that determines the asymptotics of the integral ro . we must find 
T G K"*, X G dP such that holds and t G Cx- We show the existence of such r, x in the next section. 

4.2.1. Existence of critical points. Recall that the normal cone Cp of a face F of a convex polytope Q C M"* 
is given by 

Cp = {u G M"* : {u,x) = sup(m, y), Va; e F} , (97) 

veQ 

so that Cp = for all points x G F. The purpose of this section is to prove the following general existence 
result: 

Lemma 4.3. Let Q be a convex polytope contained in E and suppose that Q (f_ dE. Then for each z G (C*)"*, 
there exists unique Tz G K™ and Xz G Q so that 

• • z) = Xz, 

• Tz G . 

Here we do not assume that Q is simple or even that Q is integral. We note that even though is 
orthogonal to dQ at Xz, the orbit {p,{e^'^^ ■ z) : r G K} is in general not orthogonal to dQ, even when Xz 
lies in a facet. (See also Figure El and the associated remark.) Our proof of Lemma lOl uses the following 
elementary, but not so well known, fact about the invertibility of Lipschitz maps. 

Lemma 4.4. |Fan | Let f : U be a Lipschitz map, where U is open in R"*. Then f has a local 

orientation-preserving Lipschitz inverse with Lipschitz constant L at a point Xq G U if and only if there 
exists £ > 0 such that 

i) liminR^o \f{x -G v) - /(a;)|/|t>| > 1/L for all x G Bs{xo), 

ii) det f'{x) > 0 for all x € Bg{xo) such that f is differentiable at x, 

iii) f{xo) ^ fidB^ixo)) and deg [/ : dB^ixo) ^ R"* \ {f{xo)}] = 1. 

Here, we let Bs{xo) = {x G R"* : |a; — xqI < denote the £-ball about xg G R™. By the degree in (iii), 
we mean the degree of ft,{dB^{xo)) G i7m-i(R'" \ {/(a^o)},^)- The hypotheses (i) and (iii) in |Fan| differ 
slightly from those above, but Fan’s proof of sufficiency uses only (i)-(iii) above. (Necessity is obvious.) 
Hypothesis (i) is given as a lemma in |h’an| . Hypothesis (iii) above is replaced in |Fan| by the condition that 
the index of / at xq is 1, which is easily seen to be equivalent to (iii) under the assumptions (i) and (ii). 

Our proof of Lemma 01 also makes use of the ‘normal bundle’ of a polytope, which we define below as 
the collection of normal cones at the points of the polytope. 

Definition: The normal bundle Af{Q) of a convex polytope Q C R™ is the subset of Tr^ = R"* x R"* 

consisting of pairs {x,v), where x G Q and v G . (Note that Af{Q) is not a fiber bundle over Q.) 

The normal bundle JV(Q) is a piecewise smooth submanifold of Tr^; it is homeomorphic to R"* via the 
‘exponential map’ 

Sq : N{Q) —> R™ , £q{x,v) = X -\- V, x G F, cG Cp (F a face of Q) . (98) 

It is easily seen that £q is a homeomorphism and is a C°° (in fact, linear) diffeomorphism on each of the 
‘pieces’ F x Cp C A/’(Q). 
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Proof of Lemma \4.S[ We note that the R™ action on (C*)"* descends via the moment map /i to an R™ 
action on S° given by r 0 fi{z) = fi{r ■ z). We let 

= M{Q) n (E° X R™) , 


and we consider the map 

We observe that the conclusion of the lemma is equivalent to the bijectivity of d), since if z G (C*)™, we have 

t ) Q x = fi{z) x = © fi{z) • z) , 

for {x,t) G M°, i.e., for x G E° and r G C§. 

To show that $ is a bijection, we let £ : E° ^ R™ be the diffeomorphism given by 

SO that 

Co fj,{w) = (log|wip,...,log|wmn • 

Thus, writing x = fJ.{w) we have 

£o$(x,r) = Q us{w) = C o [i{e^■ w) 

= (ri+log|uiip,...,rm+log|wmn 

= T + £ o ii{w) = r + £(x) . (100) 

We first show that £ o $ : M° R™ is proper: suppose on the contrary that the sequence {(x",r"’)} 
is unbounded in M°, but £ o $(x”,t") = £(x") + r” —> a G R”^. By passing to a subsequence, we can 
assume that x” ^ x° G Q- Then x° G 9E, since otherwise r" —> a — £(x°). Write = r'^u^, where 
r" > 0, \u^\ = 1. We can assume without loss of generality that u" ^ vP. We first consider the case where 
If x° = 0, then > 0 (1 < j < m), since otherwise £ o $(x",t") would diverge. But vP G C'{o}; 
hence Q C {x : {vP, a^) < 0} C R"*\E°, a contradiction. Now suppose that x? = • • • = x° = 0, x° > 0 for 
k < I < m. Then we conclude as before that > 0 for 1 < j < fc. Furthermore, in this case, rf +oo and 
r" ^ G R and thus — 0 for k < I < m, and we again obtain a contradiction. Finally, if ^ x° = 1, we 
can suppose that x® 7^ 0 and make the coordinate change xi = 1 — JP Xj , X 2 = X 2 , ■ ■ ■, Xm = Xm (which 
corresponds to permuting the homogeneous coordinates in CP™) to reduce to the previous case. Therefore, 
£ o $ : J\f° R™ is a proper map. 

Let £ = £q ■. N{Q) ^ R™, and let U = £{M°) C R™. We consider the map f : U ^ R™ given by 
/ o £\j^o = £ o $, i.e., by the commutative diagram: 


e 

U 


f 


E° 

c 

R™ 


Since £ o <I> is a proper map, / is also proper. Hence to show that d) is a bijection, it suffices to show that / 
is a local homeomorphism and is therefore a (global) homeomorphism. 

To describe the map /, for each x G E°, we let denote the point in Q that is closest to x. We note that 
Qx G E°; if X ^ Q, then Qx G dQ; if x G Q, then qx = x. Furthermore, £“^(x) = {qx, x — qx) and hence 

/(x) =x - qx+ C{qx) ■ 

We shall show that / satisfies the hypotheses of Lemma f4.41 Let F be an arbitrary face of Q n E°. To 
verify (i) and (ii), it suffices to show that detDf > 0 on the (noncompact) polyhedron 

Up ■■= £{F xCf)PU = £[{F n E°) x Cp] ■ 

To compute the determinant, we let Tp C R™, Np = Tp C R™ denote the tangent and normal spaces, 
respectively, of F. Let x° G F be fixed. For y G Up, we have 

y = x + v C{x) + V , x — x° G Tp, v G Np . 
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Choose orthonormal bases {Yr+i,... ,Ym} of Tp, Np, respectively. We let Df denote the 

matrix of the derivative {flupY respect to the basis {Yi,..., Ym} of R™. We have: 


Df = 


r‘£'(x)T 0 \ 

* I) 


where T is the m x r matrix [hi • • • W] ■ We have by Ibbll , 


( 101 ) 


1 1 

=- V 6 ^—, To = 1 - ^ Tj > 0 , (102) 

\ ^ jk Xq Xj 

for X € Ti°. Hence C {x) is a positive definite symmetric matrix, it being the sum of a semipositive matrix 
(all of whose entries are ^) and a positive definite diagonal matrix. Therefore, T^C'{x)T is positive definite, 
and hence detZl/(T) = det(T‘£'(T)T) > 0, completing the proof that hypotheses (i) and (ii) of Lemma lOI 
are satisfied. 

Note that (110211 implies that the eigenvalues of C {x) are > 1 for t € S° and hence by (IIOIII . Df{x) is a 
diagonalizable matrix whose eigenvalues are real and > 1. 

We verify (iii) by a homotopy argument: Choose a point G Q° . We contract Q to q^; i.e., for 0 < t < 1, 
we let 

Qt = {I - t){q°} + tQ , 

SO that Qo = {g°}, Qi = Q. For 0 < t < 1, Qt C E°, and hence we have a map 

: Af(Qt) ^ . 

For 0 < t < 1, we define ft : R"* ^ R™ by the commutative diagram: 

N{Qt) ^ 


R"* 


ft 


«iL 

R"* 


As before we have 

ft{x)=x-ql. + C{ql.), (103) 

where g* is the (unique) point of Qt closest to x. The above argument shows that the maps ft also satisfy 
(i) and (ii) of Ijemma Td.dl 
We write 

£ : (R™ X [0,1)) U ([/ X {1}) ^ R"* , {x,t) ^ Mx) , 

where fi=f'U^ R™. One easily sees that 



x'\ + \t-t'\ , 


and hence F is continuous. Furthermore, F is uniformly continuous on R™ x [0,t], for each f < 1. 

Let F[ denote the set of t G [0,1) such that ft : R"* ^ R"* is an orientation preserving homeomorphism. 
We shall show that iL = [0,1). We first observe that 0 G H since 

fo{x) = X — q + C{q) = T + constant . 

Note that IIIOIII also holds for the maps ft and hence as above, the eigenvalues of Dft are real and > 1. 
Hence ft satisfies conditions (i) and (ii) of Lemma [4.41 with L = 1, for all t G [0, Ij. 

We now show that F[ is open in [0,1). Suppose Q G H. Then by Lemma [4.41 Lip(/(“^) < 1 and hence 

|t - To! = 1 ^ \fto{x) - fto{xo)\ > 1 . 

By the uniform continuity of F on R"* x [0,1^^], we can choose e > 0 such that for all xq G R"*, we have 

\ftix) - ftixoY > i , for T G dBi{xo), \t-to\<e . (104) 

To simplify notation, we shall write 

deg(/,a;o,r) := deg [/ : dBr{xo) ^ R™ \ {/(tq)}] . 
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We conclude from 110411 that for \t — to\ < e, 

deg{ft,xo,l) = deg(/to,a;o,l)) = 1 , 

and hence ft is a local homeomorphism at Xq, by Lemma 14.41 By the very first part of the argument 
that f : U ^ K™ is proper, we easily see that ft is proper. Since Xq G K™ is arbitrary, it follows that 
ft : M"* ^ M"* is a covering map and therefore is a homeomorphism. 

Next we show that H is closed in [0,1) and hence H — [0,1). Let tn G H such that > to G [0,1). 
Since ftg satisfies condition (i) of Lemma lOl we can choose e > 0 so that ftgixo) ^ ftg{dBi.{xo)). Then for 
n sufficiently large, 

deg{fto,Xo,e) = deg(/t„,a;o,e) = 1 . 

It follows as above that ftg : K"* —> K"* is a homeomorphism. 

We have shown that ft is a homeomorphism for 0 < t < 1. To complete the proof of the lemma, we must 
show that / = /i is a local homeomorphism. So we let xq G U he arbitrary, and we choose £ > 0 such that 
Bi;{xo) C U and f{xo) ^ f{dBs{xf)). Then for < < 1 sufficiently close to 1, we have as before 

deg(/,a;o,e) = deg(/t, a;o, e) = 1 . 

Thus by Lemma lOl / is a local homeomorphism. This completes the proof of Lemma lOl □ 


We have thus shown (as a consequence of (PI and Lemma 14.,SI with Q =■ ip) that for each z G (C*)™, 
there exist (unique) G and g(z) = pxz G P such that the phase function 4>(-, qrz, z) : T*” x P ^ C 

given by IHHI takes its maximum real part at the point (0, q{z)), and furthermore (0, q(z)) is a critical point 
of where P is the face of P containing q{z). The latter statement follows from (P^ . (P^ and the 

fact that vectors in the normal cone Cp are orthogonal to P; i.e., Cp C F^. 

By the boundary of the normal cone Cp, we mean the boundary dCp of Cf in F-^; the complement 
Cp \ dCp is called the interior of Cp. We note the following consequence of the proof of Lemma [4. SI 

Lemma 4.5. A point z G (C*)™ is a transition point if and only if Tz G dC^f^^y 

Proof. It follows from (IIOOII that 2 Log (Tip) = P(P) + Cp, where Log (z) = (log |zi|,..., log \zm\) and that 

d[2 Log {Up)] = 2 Log {dUp) = C{F) + dCp . (105) 


Thus 


z G dTZp 4=^ 2 Log (z) G C[F) + dCp 4=^ 2 Log (z) = C{xz) +Tz, Tz G dCp . 


□ 


Remark: The identity CHS) tells us that we can decompose R"* as the disjoint union of the sets C{F) + Cp. 
If z is a transition point, then 2 Log (z) must lie in the common boundary of at least two of the sets C{F)+Cp. 
Figure El below shows the transition points in log coordinates (as solid lines) for the case where P is the 
polytope of Figure El Note that although Cpi is orthogonal to P' in Figure El Cpi is not orthogonal to 
C{F'). Similarly, in Lemma [4.51 if q[z) G F', then the orbit {/r(e’’'^^ • z) : r G M} is not orthogonal to P'. 


4.2.2. Continuation of the proof of Theorem mi Since the critical point of on the contour T™ always 
occurs on the boundary of P when z is in the forbidden region, we shall use the method of stationary phase 
for complex oscillatory integrals on quadrants. Since this general method is not so well known, we state here 
the basic result we need: 

Lemma 4.6. Let A{f,y) G C°“(]R>g x R®) such that (iy4'(0, 0) = {), A has compact support and 

2) dyd>{0,y) 7^ 0 for (0,y) G Supp (A) \ {0}, 

3) det7fj^4'(0, 0) 0 (where Tty denotes the Hessian with respect to y), 

4) < 3fi4'(0, 0) onSupp(A). 
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log 1221 

A 



l0g|2l| 


Figure 6. Transition points in log coordinates for the polytope of Figure El 


Then 


d^dy = + ciN-^ + + • • • + + O(fV-'-i)] 

for Z = 1, 2,3,..., where 


/BO -/R^o 


Co = 


(27r)'’/2yl 


\/det(-7fy«') 


«=y=o 


(If s = 0, then hypotheses (2) and (3) are vacuous, and the determinant in the formula for cq is 1.) 
Proof. (See also jSSYl-) Let us first consider the case k = 1. Integrating by parts, 




N'i’ 


A 




dy - 


«i=0 


N 



■ A ' 

Jr‘ Jo '9'Ci 



d^i dy . (106) 


Applying the stationary phase expansion jHbl Th. 7.7.5] to the first term of HI 06|l and iterating, we obtain 
the desired expansion. 

The case fc > 1 follows by induction on k, integrating by parts as above. □ 


We consider a point z in the classically forbidden region. Let Tz,Xz be as in Lemma^31 and let q{z) = pxz- 
Hence t = Tz,x = q{z) satisfies the critical point equation (li^ . and as we observed above, q{z) S dP. The 
decay function bp{z) from (I Kill is given by 

6p(2) = -vE-(0,9(2)). (107) 

Recalling lliHH) . we further decompose the integral into two parts: 

I'n = N'^ ! f e^'^^‘^’^’^^'^'>p{x)A{^,x,N;Tz)dxd^, (108) 

Jr^m Jp 

where A{ip, x, N; Tz) is given by (I92II and p is supported in a small neighborhood of q(z) and is = 1 near q{z). 
Note that it follows from (inni-iinni) that d,^d/(0, x-,Tz,z) ^ Q for x ^ q{z). Since d,p'^ does not vanish on the 
support of 1 — p{x) and sup-pm^p^' = 'I'(0, 9(2)) = —bp{z), we conclude by performing the p integration 
first, that 

X'f, = • 0 {N-°°) 

(see [Hop. Hence, we need consider only 


(109) 
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Since we are not assuming that P has interior, we let n = dim P and we let (P) denote the n-plane 
containing P. The polytope P can be given by 

P = {x £ {P) : £j{x) > 0 for 1 < j < d} , £j{x) = {uj,x) + Xj , 

where the uj are primitive vectors tangent to P. Suppose that q{z) is in a face F G P oi dimension r. After 
permuting the indices of the normal vectors {uj}, we may assume that 

F = {x G P : ij{x) = 0, for 1 < j < n — r} , 

We consider the translated cone 

r := {a: e (P) : £jix) > 0} D P 

with vertex q{z), which coincides with P in a sufficiently small neighborhood of q{z). We then choose 
(non-homogeneous) linear functions ti,... ,tr on R™ such that the map 

L = {£i,...,£ri-r,ti,...,U) : {P) M" 


is bijective and L{q{z)) = 0. We note that L : T ^ ^ L : F ^ {0} x M''. 

We shall evaluate by applying Lemma lTIHI with = £j{x) for 1 < j < n—r and y = (ti,... ... (pm) 

and with the phase 

y) ■■= (p, ti,---, tr); T, z) 

and amplitude p{x)A{ip, x, N; t). Since dy,^'(0, q{z); Tz, z) = 0 and Tz is perpendicular to P, it follows (using 
(IDlIl l that dj,d>(0,0) = 0. Since Tz is in the normal cone at q{z), the real part of the phase takes its 
maximum at q{z), and hence hypothesis (4) of the lemma holds. Also, hypothesis (2) is satisfied since 
d^'^{0,x-,Tz, z) ^ 0 for X ^ and hence dy'^{0,y) ^ 0 for y ^ 0. 

We now verify (3): Recalling 18411 . we note that the Hessian matrix (with respect to Pjx) of 4* is given 

by 

X] r, z) = H4'a(<P, X] • z) , (110) 


and hence 


where 


^^1(0,X) 


C il 

il 0 


Cjk — 


02^ 

difjdipk 


(0,x) 


p{-IjS^ + I,h) , 


(/i,...,/m) =y(e ■ z) . 


(Ill) 


( 112 ) 


Lemma 4.7. The matrix C given by nw is strictly negative definite. 


Proof. Suppose that A = (Ai,...,Am) G R™ \ {0} is arbitrary. Let v = (vi,... ,Vm), w = {wi,...,Wm) 
be given by Vj = ^/ij, Wj = ^/Yj Xj. Since ■ z G (C*)"^, Ij ^ 0 for 1 < j < m, and hence w 0. 

Furthermore ||r!|p = < 1- Therefore 

.. m 

- ‘^JkXjXk = -Ikf + (v ■ wf < -||vf llwf + (?; • r;)2 < 0 . 

^ A U — 1 


Changing variables, we see that 7fy'I'(0,0) is of the form 

0 iB 

iB* C 

where B is an r x m matrix of maximal rank r and C is nonsingular; hence Tdy 4'(0,0) is nonsingular, verifying 
hypothesis (3). 

To verify (1), we must assume that z is not a transition point. By Lemma El)! Tz is not in the boundary 
of the normal cone to F, and hence 

{tz,x — q{z)) <0 for all X e P \ P . 
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Let = L ^{E^) G F, for 1 < j < n — r (where S-j^) G R"); i.e., lies on the j-th edge of 

r. Let rjj = q^ — q{z). Then 

d^/d£,j = ^' = (r^, iij) <0 (1 < j < n - r) , 

and hence hypothesis (1) holds. 

Therefore, by (EOl), Lemma (with k = n — r, s = m + r) and nrm . 

^^^-1 + . . . + c/iV' + 0{N-^-^)] . (113) 

We let 

hp{z) = - 4 '( 0 ,g(z);Tz,z) = -{q{z),T^) +p\og 


1 + ii^r 

.1 + lie—-/2 -zP 

as in ifTHll . 

We now show that decays faster than X ]^: Recall that 
1 r /l + y”l |z,P\^^ 

= TTTvX / ^ 2 (v?) - ^7^.. . Xjvp(e^"+*‘^)dV2 , ^ 2 {p) = 0 for distx™(</?,0) > e 


(27r)^ 


1 + IPII 


where we write Tz = (ti, ..., Tm)- We choose A = A(£, z) > 0 such that 


m —ta — 


-I I '-r-^m 

1 + E,=ie 




1+ P 


< e 


l + Er=ie-^^>. 

1 + NP 


for distx'" (</5,0) > e . 


Therefore, 


Furthermore, 




\m < 


1+ u 


XNp{e^") ■ 


XNpie^^) = < #(iVP) , 


(114) 

(115) 


aGNP 


since a/N G P and is in the normal cone to P at q{z). Combining H114I) and (I115II . we obtain 

\^%\ < (iVp)”e-^P^e-^*'^(^) . (116) 

The asymptotic expansion of part (ii) of Theorem 14.1 I f for P simple) with bp(z) given by 11611 now follows 
from p?|i . fnni and dH). 


4.2.3. The decay function. In this section we verify that the function bp given by llTCll satishes parts (c) and 
(e) of Theorem I4.1f iil for any convex polytope P. 

We first note that when z G Alp, we have Tz = 0 and hence if we extend llTT)ll to all of (C*)™, we then 
have 6p = 0 on Alp. Furthermore, since z is easily seen to be continuous on (C*)'" and C°° on each 

Pp, the same holds for bp{z). We now show that bp is P by computing its derivative. Recalling that 
bp{z) = —41(0, g(z); Tz, z), we have 

-dbp = [dx4'(0, x; r, z) • Dq{z) + dz4'(0, x; r, z) • Dtz + dz4'(0, x; r, z)] . 

Here, Dtz and Dq{z) are only piecewise continuous, being discontinuous at transition points. However, from 
(EU and the fact that (0, q{z)) is a critical point of iP, we have 

dz4'(0,a;;T,z)|{2,=g(z),z=r4 = -id^'^{ip,X]T,z)\{^=Q^^=ql^z),r=r,} = 0 • 

Furthermore, since dx'^\ip=Q,T=T^ = Tz and is perpendicular to the face containing q{z), we also have 

da;4'(0, X] r, z)\{^=q(^z),r=T,} ' Dq{z) = Tz ■ Dq{z) = 0 . 

Therefore 

dbp = -44'(0,a;;r,z)|{,=,(z),z=z4 GC°((C*)'") , (117) 

completing the proof of (e). 

Property (c) is a special case of the following lemma on the dependence of the decay function 6p(z) on 
the polytope. 










34 


BERNARD SHIFFMAN AND STEVE ZELDITCH 


1 + I|Z|| 


1 + lie—W/2.^||2 


, XGE° 


(118) 


Lemma 4.8. Let P' C P C pS be convex polytopes in K™. Suppose that z G (C*)"* and let q{z) G P be given 
as in Then 

hp,{z) > bp{z) , 

with equality if and only if q{z) G P' . 

Proof. Fix z G (C*)"* and consider the function 

b{x\z) := -'E'(0,a;;T(a;),z) = -{x,t{x)) +plog 
where t{x) G K."* is given by 

.z)=x. (119) 

Note that T(q(z)) = Tz and b(q(z); z) = bp{z). Hence it suffices to show that b{x] z) > b{q{z); z) for all 
xGP\{qiz)}. 

Equation H1 1 911 yields -, and therefore 

Tj{x) = log + l0g(l -'^Xk)- \0gXj . 


( 120 ) 


Since 


dr I T—T(a:)} ■ ^!p^\{(p—Q,r—T{x)} X pp(^e 11/ ‘ '^) 0 ; 


we conclude from (tTTHll that 
and hence 


dxb — 0,T—r(a;)} — "^(^) 


d'^b 


Otu 


dxjdxk dxj 1 — X)™ 1 


( 121 ) 


( 122 ) 


Therefore, the Hessian ( gxk ) positive definite. (It coincides with the derivative C'{x) from H4.2.1I ') 
Now suppose that x° G P, x^ ^ q{z). Let 

f{t) =b{q{z)+tv,z) , v = x°-q{z). 

Then by (I121II . /'(O) = —{v,Tz) > 0, and by the positivity of the Hessian, f''{t) > 0. Therefore b{x^\z) — 
6(<7(z);z) = /(l)-/(0)>0. □ 

Now let z G (C*)™ \ Ap. We apply Lemma ^21 with P', P replaced with P, pS, respectively. Since 
qpY.{.z) = ppi{z) ^ P by assumption, we conclude by Lemma [4.81 that 

bp{z) > bpY.{z) = 0 . 

Thus we have shown that the function bp satisfies (c) and (e) for all convex polytopes P. 

4.2.4. Precise asymptotics on the classically allowed region. In this section, we prove part (i) of Theorem l4.ll 
We let P be a convex integral polytope (not necessarily simple). By ll^ . 

and we then have by 15811 . 

m 

n|Arp(2,z) = Y[{Np + j) - Rn{z) , 

where 


1=1 


1=1 


1 


IIXo 


rlx^"(^)P • 


Rn{z) = 

aGNpS\NP 

To simplify our computations, we introduce the renormalized monomials 


~Np 


(z) := 


(iVp)! 


_{Np + my.\ llx. 


Xa{z) _ (Np 
a 




Np 

a J (l + ||zP)^p/2 


(123) 


<Np. (124) 
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Let a G NpYj \ NP be fixed. We easily check that 

z) = ^ 


{Np + to)! 


{2'kY‘ 


g'I'iv(<p,a;'r,z) 


d(p , 


(125) 




(126) 


where the phase is given by 

(Equations 1125II " (I126II also follow from (El and m with P = {a} and p replaced by Np\ 'i’N is the phase 
with p replaced by Np.) Since 3(?'l'Ar((p, x; r, z) < 'I'Ar(0, x; t , z), we have 

|m^^(2)P < . (127) 

Choosing r such that Npp{e~'^^^ ■ z) = a and recalling (I118II . we have 

-«'Ar(0,a;r,z) =7V6{-}(z) . (128) 

By Lemma ITHl > bp{z) = 0 for a; € pS \ P° and z € Ap (i.e. pp{z) G P°). 

Now let iL be a compact subset of ^p, and let 

Aa = i inf{6{a;}(z) :xGpJ:\P°, z G K} > 0 . 

Therefore, 

|m^^(z)p < , for z G iL, a G NpY. \NP, N> 1. 

Since Y^iNp'S \ NP) < (Np)"^, we conclude from 112811 and 112911 that ||i?Ar||c“(A') = 0{e~^^^). 

To verify the estimate, we recall from EU that 

'Ap/.M2 /'Np') 


hiz)^\zY^ , h{z) = (1 + ||zf )-P . 
Differentiating with respect to the variables pj = log \ zj\, we obtain 

dpi\m^Piz)\^) = (^^)hiz)^\zY^[2a + Ndplogh] = |TO^^’(z)|y2a + iVdplog/i] , 


(129) 


and hence the estimate follows from the C° estimate. Differentiating repeatedly, we obtain all the 
estimates, completing the proof of part (i) of Theorem 14.11 


4.2.5. Non-simple polytopes. In this section, we prove the general case of part (ii) of Theorem 14.11 by a 
reduction to the simple case. 

Let P be a non-simple convex integral polytope of dimension n. Consider a face P of P of dimension r 
and let z G TZp. We say that P is simple at P if //J{F) = n — r. (Recall the notation from H2.ll Note that 
we always have //J' (P) > n — r. A polytope is simple <;=^ it is simple at each of its vertices it is 

simple at each of its faces.) 

We first consider the case where P is simple at P. (This is always the case when codim P < 2.) We 
construct a simple integral polytope Q D P coinciding with P near q{z) as follows: Recall that q{z) G P 
(by the definition of TZf), and consider the ‘barrier cone’ Bp/P) generated by elements of the form x — q(z), 
where x G P. Choose L G GL(to,Q) such that 

L : Bf(P) « R”-’’ X R’' X {O}™"” . 

(E.g., let L map the interior normals {uj : j G Jf} at F to the first n — r standard basis vectors in R™ and 
map a basis for the tangent space of F to the next r standard basis vectors.) We then let 

Q = q{z) + L-^{[0, MY-'- X [-M, MY X {O}™"’^) , 

where M G is chosen so that the vertices of Q lie in . By the construction, Bf{Q) = Bf{P) and 
hence Q coincides with P near q{z). We then replace M with a sufficiently high multiple if necessary, so 
that Q G) P. 
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Since the integral polytope Q may not be contained in pS, we consider the convex rational polytope 
Q' := Q npE. Recalling 1124II . we have 

^Np{z,z) = nij^Q,{z,z) - ’ 

where by the Szego kernel IIi^q/, we mean where Kn is the convex hull of HNQ'. By 1127I1 - (I128II . 

for a € {NQ' \ NP) n Z"*, we have 

< exp {-Nb[a/N}iz)) < , 

where 

C = inf{5{a;}(z) :x&Q'\P}. 

Since Q' coincides with P at q{z)^ Tz is in the normal cone to Q' at q{z) and hence q{z) = qQ'{z). Therefore 
by Lemma 14.81 b{x}{z ) > bQ'{z) = bp{z) for all x G ^ q{z), and hence C > bp{z). Since #{NQ') < 

it follows that 

\\^\Np{z,z) — W-\p[Qi{z,z)\\co(^K) = 0{e bp{z) < C'< C, K az TZ°p. 

The similar C* estimates follow as in the proof of Theorem 14.iT il in 

Thus, to obtain the asymptotic expansion of Theorem l4.ir iil for P, it suffices to obtain the same expansion 
for Q'. To do this, we make the simple observation that 

^NQ'iz,z) = ■ z)xNQ'{e"^)dip = ■ ^)AwQ(e*^) , (130) 

since lattice points in NQ \ NpT, = NQ\NQ' do not contribute to the integral. By repeating word for word 
the argument at the beginning of ^4.21 and in H4.2.2I with (118011 in place of II70II , we obtain the expansion 

U\NQ'iz, z) = K(^) + cf {z)N-^ + • • • + cf {z)N-'^ + (z)], 

with Cq > 0, and hence part (ii) of Theorem 14. II holds for P. 

Now suppose that P is not simple at F. We subdivide the barrier cone Bp{Q) into simple rational cones 
: fc = 1,..., s} so that D Bp{F) for all k and intersections of any number of the B^ are faces of each 
of them. (Note that Bp{F) = Tp.) For a non-empty subset £ of {1,..., s}, let Pg = P D Hfeef (9('^) 
which is simple at F. Since Tz is in the interior of the normal cone of P at q{z), Tz is also in the interior 
of the normal cone at q{z) of each of the Ps- (In particular, z is not a transition point for any of the Ps-) 
Thus by the case proven above (which does not use the fact that the vertices of P are lattice points, only 
that the vertices of Q' lie in Z"*), we have the expansion 

n\^P^{z,z) = + 4iz)N-^ + ■ ■ ■ + cUz)N-^ + R^iz)] , z G 7^^ , (131) 

where Cq > 0. Note that 6p(z) is the same for all of the Pg. Then by the inclusion-exclusion principle 

n|^p(z,z) = ^(-l)#^-in|^p,(z,z) (0^£c{l,...,s}). (132) 

£ 

The desired asymptotic expansion for n|^p(z, z) now follows from (I131II - (I132II . by noting that n|^p(z, z) > 

IIijvp^ij (z, z) and hence c^(z) > cl^\z) > 0. 

The proof of Theorem EU is now complete. 


4.2.6. Proof of Proposition \4 - Bv Theorem l4.lT iil. 

^logn|Arp(z,z) ^ -bp{z) , 

for all non-transition points z. We use the log coordinates pj + iOj = logZj, so that 

u^(z) = llogn|^p(z,z)+plog(l + ||zf) = llog ^ zG(C*) 

a£NP 


Thus for all non-transition points z, we have 

un{z) Uoo{z) := plog(l -b ||z|p) - bp{z) . 


(133) 

(134) 


( 135 ) 
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We must show that convergence of (insj also holds at the transition points and is uniform on compact 
subsets of (C*)™. From (Tmil . we obtain 

J .. . J lNp\„2(a,p) l^aeNP \ a ^ 

^dplog 2^ („)e ' - (Np) 2 {c.,p) ■ 

aeNP l^adNPy a 

Since ||a|| < Np for all a G 7VP, we therefore have the uniform upper bound 


(136) 


II pWa||£”((C*)'") — fJY' ('^P)p2(a,p) — 

l^aeNP I a 


Since {mjv(z*^)} converges for any non-transition point it follows that {mat} is uniformly bounded and 
uniformly equicontinuous on compact sets. Therefore it converges uniformly on compact sets in (C*)"*. □ 


5. Distribution of zeros 


5.1. Expected zero current and the Szegd kernel. In order to deduce Theorems ll.3tn~^ from Proposi¬ 
tion oi we need to relate the expected zero current to the conditional Szego kernel n|p. This relationship 
is given by the following result: 

Proposition 5 . 1 . The expected zero current ofk independent random polynomials fj G Poly(Pj), 1 < J < fc, 
is given by 

k k 

E|P„...,pJ%....,/J = A E|P,(^/,) = A 

i=i i=i 

Recall that wps = ■^'^^591og(l -I- ||7:||)^. The case fc = 1 of the proposition is essentially the same as in 
ISZII Prop. 3.1], but neither the statement nor the proof there cover the application we need. We remark 
that Proposition IQ is a special case of a more general statement for the zeros of general linear systems on 
compact Kahler manifolds, which we discuss in jSZ.3j . 

We recall that the space of (fc, ^)-currents on an n-dimensional complex manifold Y consists of 

the continuous linear functionals on the space of C°° compactly-supported {n — k,n — Z)-forms 

on Y. We say that a sequence j G V'^'^iY) converges weakly to ^<0 if {'^o,(p) for all test forms 

ip G 

We begin the proof with the following lemma, which covers the k = 1 case. 




V 27r 


ddlogIl\p.{z,z)+pjU;Fs] on (C*) 


Lemma 5 . 2 . Let P be a convex polytope in pY (or more generally, an arbitrary subset of pY n Z"*/ Let Y 
be an algebraic submanifold o/(C*)'", and let h G C°°(Y) be given by 

h{z) = logn|p(z,z) -f plog(l -f ||zf) , 

Then 


z G y. 


E|p(-^/|y) — 


271 


-ddh G V^’^iY) . 


Proof. The proof is essentially the same as the proof of Proposition 3.1 in |SZ1| . Write fy = f\Y for 
/ G Poly(P). li fy ^ 0, then the current of integration over the zeros of / is given by the Poincare-Lelong 
formula: _ 

Zf^ = ^dd\og\fy \G V'^’\Y) . (137) 


Recalling m , we can write 


/- = E 


aGP 


IIXo 


rXa = (c, G) 


' = (cc) 


aeP 


G = 


llxo 


rXalY 


aGP 


(The norms ||xa|| are computed over all of C™ as before; see (imi - (ll7)ll .l We then write G{z) = |G(z)|u(z) 
so that |u| = 1 and 


log I (c,G) I = log|G| -blog|(c,u)| . 


(138) 
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By (II ,S7ll - n.SSII . we have for any test form ip G P" ^(F) (where n = dimy), 


lE\p{ZfY),ip) 




'C'' 


log|(c, G)| d'y{c),dd(p 


dd / log|G|(i7(c),(/j 
Jc’‘ 




dd / log|(c,M)|d7(c), ip 


\ Jc^ 


where k = and d'y{c) = -^e dc. Upon integration in c, the second term becomes constant in z and 
dd kills it. The first term is independent of c so we may remove the Gaussian integral. Thus 


E|p(Z/,) = ^aaiog|G|^ 


(139) 


Recalling lITTI) and (IHHi . we have 

n|p(z,z) = (i + ||zf)-^'|G(z)p , zer, 

and the formula of the lemma follows. □ 


Proof of Provosition Um We first note that the expected current is well dehned, since for almost all choices 
of the fj, 

[ (^FS ’'< Pi ■■■ Pk , 

and hence for each test form p G x>rn-k,m-k function (/i,..., fk) , v?) is in (Poly(Pi) x 

• • • X Poly(Pfc)). 

We shall verify the current identity by induction on k. Proposition 15.21 with Y = (C*)™ gives the case 
k = 1, so assume that fc > 1 and the proposition has been verified for fc — 1 polynomials. By the inductive 
assumption, it suffices to show that 

= E|p^_..._p,,_,^(Z/^,... A E|pJZ/J . (140) 

Write Y = |- By Bertini’s Theorem, we know that the Zf. are smooth and intersect transver- 

sally in (C*)™ and hence Y is smooth (and of codimension k— 1) in (C*)™ for almost all /i,..., fk- Therefore 
for p G .^^g 


iZh,...,h,p) = <p= p = {Zf^\Y,P\Y) 

JZf, f. dYnZf^ 


.A 


(141) 


for almost all fi,..., fk- Averaging II141II over fk and applying Lemma [5.21 we obtain 

/ (%,..../fc)'P)^7pfc|Pfc(/fc) = E|pjy/j^|y,(^|v) = / E|pjZ/J A (/J = E|pjy/J A V?) , 

JPolv(Pt) Jy 

(142) 


/Poly(Pfc) 

where 


27T 


E|P.(%) = ^logn|p,(z,z)+pu;FS G2?'’'((C*)™) . 


Averaging 114111 over /i,..., fk-i and applying 114211 . we obtain 

Jk^P) = '^\Pi,...,Pk-i{Zh,...jk-i, E|Pj^(Z/JA(/j) = (E|Pj_...^pj^_^(Z/i,... jj,_jAE|pjZ/J, p) . 

□ 
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5.2. Asymptotics of zeros: proofs of Theorems O d.nd To prove Tlieorem^^^^ avg let Uoq 
be as in the proof of Proposition 14.21 By Proposition 15.11 (recalling that wfs = ■^^^591og(l + ||2:|P)), we 
have 

y^dduN = ^^NpiZf) on (er, (143) 

and hence un is plurisubharmonic. By ini and dSl), we conclude that 

^E|jvp(Z/) ^ = PUJFS - ’ (144) 

where differentiation is in the sense of currents (see e.g., |Shl Chapter I]). 

We now show that the current ddbp G is given by a (1, l)-form with piecewise smooth 

coefficients; in fact ddbp is C°° on each of the regions TZp given by fl^ . (Equivalently, the Radon measure 
ddbp A tups does not charge the set of transition points, and hence formula (tnii can be interpreted as 
differentiation in the ordinary sense on the regions TZp.) By Theorem 14.11 bp G C^((C*)'"); i.e., if z® G 
dTZp n dTZpi, then the values of dbp{z^) computed in the two regions TZp and TZp' agree. Then for a test 
form ip G 2?™-b"^-i((C*)™), we have 

{ddbp, p) = ^ f bpddp = ^ f ddbp Ap — ^ f (dbp Ap + bp A dp) . (145) 

p' J 7^ p F " ^ p F dF- F 


Note that dTZp consists of C°° real hypersurfaces (consisting of those points of dTZp that are contained in 
the boundary of only one other region TZpi) together with submanifolds of real codimension > 2, and hence 
Stokes’ Theorem applies (see e.g. 4.2.14]). Since bp and dbp are continuous on (C*)™, the boundary 
integral terms in 114511 cancel out and we obtain 


{ddbp, p) = ( ddbp A p = f ddbp A p , 

J'R% J{C-)^\E 


(146) 


where E := IJ^ dTZp is the set of transition points. Formula 1 14611 says that the current ddbp is a (1, l)-form 
with piecewise smooth coefficients obtained by differentiating bp on the regions TZp. 

We now let 


ipP = dduoa = pwps — -^^—ddbp 


(147) 


271 - ^ 27r 

on each of the regions TZp. By (jnni), the current ij^p is also a piecewise smooth (1, l)-form; by 114411 

A“^E|jvp(^/) ^ i/'p (weakly) . (148) 

We shall show convergence of (ITlHll when we prove Theorem n~i below. Continuing with the proof 
of Theorem o we observe that (ii) is an immediate consequence of (tTTTIl . since 6p = 0 on the classically 
allowed region. 

To verify (iii), we again use the log coordinates Q = pj + iOj = logZj, so that 


Moo = plog (^1 + ^ - bp{eP) . 

Since Uao is independent of the angle variables 6j, we have 

Ipp = ^^^^^dduao = X! ^0 AdC,k>0 . 




87 ^ ^ ^Podpk 


(149) 


(150) 


Thus we must show that the Hessian of Uoo{p) has rank r at points G TZ°p, where r = dimF. From 11 1 711 
and inu, we have 


dpbp = -pdplog (l + ^ 

To simplify (tTCTl . we note that for a G K™ we have 

2g-<Ti+2pi 


dplog(l + ^e 


-pdplog(^l + ^ 

2g-fTi+2pi 


= 2P3 


1 ^ g-CTj+2pj ’ ■ ■ ■ ’ X ^ g-<Tj+2p3 


(151) 

= 2p{e-°l'^ ■ z) . (152) 
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Hence by (tTK^ with a = Tz and 0, we have 

dpbp = 2p{fi - , 

where ^ = (^i,... ,^m) ■ (C*)"* ^ (C*)"* denotes the map given by 
By (1133, we then obtain 

dpUoo = 2pfio^. 

Hence the Hessian TipUoo equals the Jacobian of 2pp o Since p o ^{z) = ^q{z), we easily see that 


poi-.n°p^-F 


(153) 

(154) 

(155) 

(156) 


so rank_D(/4o^) < r. In fact, II156II is a submersion, so that the rank equals r. To see that it is a submersion, 
we recall from the proof of Lemma lT!^ that the ‘lipeomorphism’ 4)“^ : E° ^ Af° restricts to a diffeomorphism 

E° D p{'R-°f) ^ X Cf C J\f° , p{z) ^ (^^q{z), Tz^ = (m o C(z), Tz) , 

and therefore (IT?^ is a submersion, completing the proof of (iii). 

We now prove Theorem 11.41 the case A: = 1 of the theorem will then yield part (i) of Theorem 11.31 Let 
Pi,..., Pfe be convex integral polytopes, as in Theorem ll.4l We first show that 

A^"''E|Arp(%,...,/J ^ V’fe := '0Pi A • • • A ipp^ weakly. 

For 1 < j < fc, we let 


(157) 


Un = logn|Arp^.(z,z) +pjl0g(l + ||zf) , 


SO that, recalling 1123, 

By (inzj, 


u%{z) ui^(z) :=Pj\og{l + \\zf) - bp^{z) . 


^ddui, = V;p, . 


Recalling Proposition 15.II we introduce the (fc, A:)-forms 


kn 


:= (A • • • A f ^dd’ 


V 27r 


\ 27T 


=iV-'=E|ivP,.....AP, (%,...,/,). 


Since locally uniformly, it follows from the Bedford-Taylor Theorem mi HD that 


Kn 




A • • • A 


\ 2tt 


ddu'l 


weakly. 


(158) 


(159) 


In fact, the limit current in (I159II is absolutely continuous, and hence is equal to the locally bounded (piecewise 
smooth) (fc, fc)-form ipk- To see this, we note that by the Bedford-Taylor Theorem, 

k ^ r /— T 1 ^ 

/\ (^P, *iPe)= /\ -^ddiui^ *ip^) -> /\ 

i=i 1=1 ^ 1=1 


2tt 


where denotes an approximate identity. Since the currents ipp^ have coefficients in the forms 
Aj=i('*Ap * Pe) have locally uniformly bounded coefficients; absolute continuity of the limit current follows. 
The weak limit (1123 now follows from n^-(tT^. 

To complete the proof of Theorem ll.4l we must show that kn —> V'fe in for all compact K C (C*)™. 

Let £ > 0 be arbitrary, and choose a nonnegative function rj G P((C*)'") such that 77 = 1 on a neighborhood 
of {El U • • • U Ek) (A K, where Ej is the set of transition points for Pj, and 


[ /\ rjujps ^ < £ ■ 
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By (I157II - 1158II . 

I KN for TV ^ 0 . 

J(C*)- 

(Note that the above integrands are nonnegative.) Therefore, 

||ka - AWc^k) < 11(1 - v){kn - i^k)\\ci{K) + \\v>^n\\ci(k) + WvMchk) 

Since is a positive (fc, A:)-form, 

IIlyKivII£1 = J K]\f A lyt^FS ^ ■ 

Similarly, ||?7i/'fe||£i(_ft:) = Jj^tpk /\ Si^^ce 

{I — r]){KN — i^k) ^ 0 uniformly on K , 

it follows that 

||kjV — '^k\\ci(^x) “*■ 1* ’ 

completing the proof of Theorems ^31 and O 


□ 


We now derive our integral formula from 115811 : 

Proposition 5.3. The decay function bp of Theorem \l.^ is given by 

bp{z) = r f-(;(e-"/2 . z) + • z) 

Jo L 


• da . 


Proof. Fix a point z in the classically forbidden region. Then by 1 15311 with the change of variables pj = 
—(7^/2 + log \zj\, we have 

bp{e~'^^^‘^ ■ z)) — bp(z) = f dcbp{e~'^^^ ■ z) ■ da = p f [p o ■ z) — ■ z)] ■ da , (160) 

Jo Jo 

where the map ^ is given by itTCll . By definition, q(z) =p^o^(z), and hence p^o^(e • z) = g(e ^I’^-z). 
Furthermore 6p(e“'^*/^ • z) = 0 since • z = ^(z) G dAp, and formula I38II then follows from (tlM . □ 


We note that the proof of Theorem n~i also gives an asymptotic expansion away from transition points: 

Theorem 5.4. Let Pi,..., Pk be convex integral polytopes. Let U be a relatively compact domain in (C*)™ 
such that U does not contain transition points for any of the Pj. Then we have a complete asymptotic 
expansion of the form 

1 , , </5i 

^E|jvPi,...,APfc(-^/i....,/J ~'0Pi A • • • A + — H-- 

with uniform C°° remainder estimates, where the pj are smooth {k,k)-forms on U. 

Proof. By 11 3411 and we have 


^ _ _ 

E|APi.....WPfc(%....,/J = /\ -99 


ATfc 


1=1 


27r 


N 


logBiATp^ (z, z) + Pj log(l + ||zf) 


The conclusion now follows from the asymptotic expansion of Theorem 14.11 


□ 
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5.3. Proof of Thfiorem 11 .1 ( Theorem o is an easy consequence of our prior results. Indeed, by Propo¬ 
sition EH we can write N ™E|7 vp(^/i,= Gnuj'^, where the Gn are positive C°° functions on (C*)™. 
We must show the convergence of the sequence of measures 

dXN := iV-™Epvp(%,..../„) = Gn^^ . 

Now let i? be a Borel subset of (C*)™. Then 

Xn{B) = iV-™E|^p(#{z G B : /i(z) = ■ • • = U{z) = 0}) . 

We first consider the case where B is contained in a compact set K C (C*)™. Then by Theorem 11.41 with 
k = m and Pi = ■ ■ ■ = Pk = P, it follows that 

Gnuj'^^iPp in C\K) . 

By Theorem 11.31 on and ipp = 0 on the complement of ^p since its rank is less than m 

there. Therefore, 

Xn{B)= [ xbGnco"^^ f xsnApP’^^FS =^!p™VolcP-(Bn^p) . (161) 

Jk J K 

For the general case, let e > 0 be arbitrary and choose B' C B such that B' is relatively compact in 
(C*)™ and Volcprr.(i3 \ B') < e. Then by the above 

liminf Aat)!?) > liminf AAr(i30 = rn\p'^\o\cv^{B' n Alp) > m!p"* Volcp™(i? H Alp) — m\p^e . 

N —>00 N —>oo 

Since e > 0 is arbitrary, we conclude that 

liminf AAr(i?) > mlp™ Volcp™(i? H Alp) . (162) 

N—>oo 

To obtain the reverse inequality, we recall from the Bernstein-Kouchnirenko Theorem EiESIlEia that 
the number of common zeros of {/i,..., fm} equals m!7V’"Vol(P) for almost all (/i,..., fm) G Poly(P)™, 
and hence 

AAr((C*)™) = m!Vol(P) for all iV > 1 . (163) 

It is well-known and easy to verify that Volcp™(/r“^(n)) = Vol(fl) for any C E. (This is a special case of 
a general fact about the moment map in symplectic geometry; see (116611 in the Appendix.) Recalling that 
Alp = ^“^(1-P°), we then have 

Vol(Alp) = ^Vol(P) . (164) 

It follows from II162II with B replaced by its complement P'^ = (C*)”^ \ B and (I163II - II164II that 

limsupAAr(P) = m!Vol(P) — liminf Aiv(P‘^) < m!p™ Volcp"*(P H Alp) . (165) 

N^oo N^oo 

Therefore, (116111 holds for a general Borel set P C (C*)™. □ 

5.4. Vanishing of ipp along the normal flow. The proof of Theorem ll. 31 1111 gives us some more informa¬ 

tion about the expected zero current in the classically forbidden region. Roughly speaking, let / G Poly(iVP) 
be a random polynomial with Newton polytope NP, with N large, and let G \Zf \ \Alp; then \Zf \ is highly 
likely to be close to being tangent to the orbit of the normal flow at In particular, if is in the flow-out 

of an edge (dimension 1 face) of P, then T^o{\Zf\) is likely to be a good approximation to the tangent space 

of the normal flow through 

To make this statement precise, we define the complexified normal cone of a face P, 

Gp := {t + iO : T G Gp, 9 G Tp} . 

(Recall that Gp C Tp.) We note that Gp is a semi-group, which acts on TZp by the rule r]{z) = e^ ■ z] 
we call this action the ‘(joint) normal flow.’ The (maximal) orbits of the normal flow are of the form 

Gp ■ z^ = {e^ ■ z^ : r] £ Cp}, where z^ G We note that the orbit Gp • is a complex (m — r)- 

dimensional submanifold (with boundary) of (C*)"*. (Indeed, (Cp • z°) (A TZp is a submanifold without 
boundary in TZ°p.) 
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Theorem 5.5. Let P be a convex integral polytope and let ipp be the limit expected zero current of Theorem 
O Then ipp vanishes along the orbits of the normal flow. 

Proof. Let 

0:=Cf-z° = : r G Cf, 6» G T^} 

be a maximal orbit of the normal flow, where pp{z^) G F. For z = • z° G O, we have 

M O az) = K{z) = -qie^ ■ Z°) = -q{z°) = p{z^) 

p p p 

and hence fj. o ^ is constant on O. It then follows from H15()ll and (I155II that ijjplo = 0. □ 

To relate Theorem 15.51 to the explanation above, suppose for example that z^ is in the flow-out of a 
1-dimensional face (edge). Choose coordinates wi,... ,Wm, so that the normal flow near z° is given by 
wi = constant. Since tpp vanishes along the orbits of the flow and has rank 1 near we have 

—'E\p[p{Zf) ^fp = ic{w)dwi A dwi , c{w) > 0, (166) 

near z°. For a regular point z of \Zf\, we let rif{z) G T*i>0((C*)'") be a unit vector (unique up to the 
action) annihilating r^’°(|Z^|). We can then write Zf = Szfi^r] /\ fj) , where Szf is the measure given by 
{5zf,<p) = J\Zf \ ‘Pdyol 2 m- 2 - Writing rj = J^jLi aj{w)dwj, we see by (11661) that 

4:E|Afp / |aj|^dVol 2 m -2 ^ 0, for j>2. 

^ J\Zf\ 

Thus the expected value of the average distance between the tangent spaces of \Zf\ and of the normal flow 
approaches 0 as iV ^ oo. 

5.5. Amoebas in the plane. The term ‘amoeba’ was introduced by Gelfand, Kapranov and Zelevinsky 
[SEZ] to refer to the image under the moment map of a zero set of polynomials, and have been 

studied in various contexts (see |FPTI ICKZI IMil I IPR| and the references in the survey article by Mikhalkin 
|Mi2| l. The image of a zero set under the moment map p is called a compact amoeba, while the image under 
the map 

Log : (C*)™ ^ M™, (zi,...,Zm) 1-^ (log|zi|,...,log|z„^|) (167) 

is a noncompact amoeba, or simply an amoeba. Note that Log is the moment map for the T"* action with 
respect to the Euclidean symplectic form '^dxj A dyj, and Log = o p,, where £ : E° « is the 
diffeomorphism given by (P|) . 

To illustrate what our statistical results can say about amoebas, we consider zero sets in (C*)^. An 
amoeba in is the image of a plane algebraic curve under the Euclidean moment map Log. An example 
of an amoeba of the form Log(Z^), where / is a quartic polynomial in two variables with (full) Newton 
polytope 4S, is given in the illustration from |Th| reproduced in Figure |7| below. 

One notices that this amoeba contains 12 ‘tentacles’. By definition, a tentacle on a compact amoeba A 
is a connected component of a small neighborhood in A of A n 9E; the tentacles of a noncompact amoeba 
correspond to those of the compact amoeba under the diffeomorphism E° R"*. 

It is known that there is a natural injective map from the set of connected components (which are convex 
sets) of the complement of a noncompact amoeba A to the set P 0 7? of lattice points of the polytope, and 
that there are amoebas for which each lattice point is assigned to a component of the complement. (This fact 
is also valid in higher dimensions; see For a generic 2-dimensional noncompact amoeba with 

polytope P, each lattice point in dP corresponds to a distinct unbounded component of R^ \ A, and adjacent 
lattice points correspond to adjacent unbounded components. (The correspondence is given in [Mil I §3.1] 
or |FPT| .l Each tentacle thus corresponds to a segment connecting 2 adjacent lattice points on dP. Hence 
the number of tentacles of A equals the number of points of dP (A Z; this number is called the length of dP. 
(For example, the 12 tentacles in Figured correspond to the 12 segments connecting the 12 lattice points in 

dm.) 

We can decompose dP into two pieces: d°P = dP Op'Z° and d'^P = P O d{pY,). Tentacles corresponding 
to segments of d°P end (in the compact picture E) at a vertex of E, and tentacles corresponding to segments 
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Figure 7. An amoeba with polytope 4S 

of d^P are free to end anywhere on the face of S containing the segment. We call the latter free tentacles, 
and we say that a free tentacle is a classically allowed tentacle if its end is in the classically allowed region 
Ap. For an amoeba A, we let vat (A) denote the number of classically allowed tentacles of A. It is clear 
from the above that 

uat(A) < #{free tentacles} = Length(5®P) 

and that this bound can be attained for any polytope P. Here, ‘Length’ means the length in the above sense; 
i.e., the diagonal face of pE is scaled to have length p. As a consequence of Theorem I l.dl (for m = 1), we 
conclude that this maximum is asymptotically the average: 

Corollary 5.6. For a convex integral polytope P in we have 

^E|Arp(uAT(Log(Z/)) ) ^ Length(a'’P) . 

Proof. Let Fi,F 2 ,Fz denote the facets of pE and apply Theorem 11.51 111 to the 1-dimensional polytopes 
PnF„j = 1,2,3. □ 


6. Appendix: approach using toric geometry 

Although our results and proofs do not depend on the theory of toric varieties, the subject motivated 
some of our ideas. We describe briefly in this appendix how toric varieties can be used to give a geometric 
derivation of the Szego kernel asymptotics of Theorem 14.11 

6.1. Relationship to toric geometry. A toric variety is a complex algebraic variety M containing the 
complex torus (C*)™ as a Zariski-dense open set such that the group action of (C*)'" extends to a (C*)"* 
action on M. A toric variety can be constructed from a fan by gluing together the affine varieties arising 
from the cones in the fan (see |Ful Chapter 1]). The toric variety Mp constructed in this way from the fan 
of a convex integral polytope P is a projective variety. If P is simple, then Mp has orbifold singularities; if 
P is Delzant, then Mp is smooth ini (see also §2.1-2.2]). 

For example, CP™ is the toric variety corresponding to the simplex E C M™. Hence, the toric variety 
corresponding to the square from Example 1 in o is the product CP^ x CP^ of projective lines. When P 
is the trapezoid in Example 2, the associated toric variety Mp is the blow up of (0,0,1) G CP^, i.e., Mp is 
the Hirzebruch surface Pi. The toric variety corresponding to the polytope in Example 3 is the Hirzebruch 
surface P„ (see |Ful §1.1]). 

We shall now assume for simplicity that P is Delzant and has nonempty interior. In this case, Mp can be 
given as the closure of the image of a monomial embedding : (C*)™ ^ (see IHkTI Chapter 5] 
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or [HtzII). We also give Mp the structure of a symplectic or Kahler manifold with symplectic/Kahler 
form Up := ihpWFS- (The symplectic form up depends on the choice of constants defining the monomial 
embedding see |STZl| .l We define the line bundle Lp <i)pO(l), where 0{\) denotes the hyperplane 
section bundle on and we give Lp the Hermitian metric obtained by pulling back the Fubini-Study 

metric on 0(1) so that Lp has curvature form up. 

For example, for the case P = S, the toric variety Ms = CP™, Ls = 0(1), and 

Poly(iVS) ^ iJ°(CP™,0(iV)) = M°(Me,L^) . 

(Recall that i/°(M, L) denotes the space of holomorphic sections of L.) More generally, 

Poly(7VP) ~ P°(Mp, L^) = . (168) 

The underlying real torus T™-action on the symplectic manifold {Mp,up) is Hamiltonian with moment 
map /rp : Mp P fibering Mp as a singular torus bundle over the polytope P (see e.g., jFul §4.2]). 
The volume form on Mp can be written in terms of the moment map: ;^Wp = dl A d9, where dl = 
fj,*{dxi A • • • dxm)- Integrating, we obtain 

Vol(P)=/ dlAd0= [ = ^deg[ci(Pp)™] , (169) 

Jmp Jmp m\ m\ 

which yields Kouchnirenko’s Theorem, since deg[ci(Pp)™] equals the number of points in the intersection 
of the divisors of m generic sections in H^{Mp, Lp) ~ Poly(P). (This proof of Kouchnirenko’s Theorem is 
from Atiyah isi.) 

Remark: Using the isomorphism CSHl), one can view xp(e*®) as the character of the torus T™ action on the 
space P°(Mp, Lp) of holomorphic sections of the line bundle Lp over the toric variety Mp (or equivalently 
as the equivariant index of the 8 operator on sections of Lp). The characters {xNp{e^^)} are the characters 
of the powers Lp of Lp and are given by the equivariant Riemann-Roch formula |BV3I l(lu| (see also EH 
§8]). When 0 = 0, the Riemann-Roch formula gives another form of the Ehrhart formula for the number of 
lattice points of NP: 

XNp{f) = #{NP) = dim Poly(7VP) = dimP°(Mp,L^) 

= x{Mp,L^) = ^deg[ci(Pp)'=UTodd™_fc(Mp)] — 
k^O 

= jv" + ■ ■ ■ + degTodd^(Mp). (170) 

777! 

(Note that equations it^ and 117011 provide an alternate derivation of 116911 . 1 


6.2. Geometric proof of Theorem l4.1L We outline here an alternative approach to evaluating the integral 
(IHH by lifting it to Mp, obtaining a complex oscillatory integral over T™ x Mp. By (IMIl . we have 


-^N — 


iV” 


( 271 )^^ 


d>p ■ Ci(:^)e^'^“(^’^’^)Todd(P, N-^d/dh) 


^N{T+iip,x) 


lp(h) 


(171) 


h=0 


We can analyze this integral by lifting to the toric variety {Mh,Uh) associated to the deformed polytope 
P(/i) given by (IHHI) . The underlying complex variety M/, = Mp is fixed, while the symplectic form is given 

by 


Uh = Up hkCi{Lk) , 

k=l 

which is affine in h. Here, L^ is the line bundle associated to the divisor fiJ,^{Fk) in Mp, and Ci(Pfc) is a 
T™-invariant (1, l)-form in the Chern class of L^. 

The moment map Hh = (if,... ,/„) : P{h) is also an affine function in h. Indeed, let Xj denote 

the vector field on Mp giving the infinitesimal action on T™ corresponding to the j-th unit vector in the 
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Lie algebra ~ K™. Then we have: 

d d 

dij = dfih{Xj) = iVhiXj, •) = ujp{Xj, *) + X! ^kCi{Lk){Xj, •) = dij + ^ hkCi{Lk){Xjp) 

k^l 

and hence 

If = I,. 

k^l 

Using the fact that 

[ f{x)dx=f 
Jp(h) JMp 


we then lift nTTl) to Mh to obtain: 


-Ln — 




Todd{T,N-^d/dh)\h=o 


I Mp 




(172) 


i+Er= 




J=1' 


(27r)"*TO! 

The phase is now the function on T”^ x Mp given by 

d 

= {r -\-i(f,fip{w)) hk{r -\-i(f, J^{w)) -\-p\og 

k^l \ 

The amplitude is a polynomial in h: 

= Wp + ^ /l“7a = Giw, h)uj^ . 

l<\a\<m 

Interchanging the order of integration and Todd differentiation in (ITTa and then using as before, we 
obtain: 


(173) 


(174) 


7"i _ 


TV" 


(27r)'"m! 


ip, w, r) u;^{w) dp 


(175) 


/T^ JMp 

with amplitude 

A{N,p,w;t) = ^i{p)G{w,N-'^d/dq)Todd{P,q)\g^^^^_^_i,^ji.(^^-^y 
The phase is our familiar phase function dt given by llilTll lifted to Mp. 

We need to find the critical points where the phase has maximal real part. As we determined before, 
the phase has maximal real part where = 0 and r is in the normal cone to P at fj,p(w). As before, we find 
that the critical point equation dip'i’{p, p,p(w);T, z)\,p-Q = 0 is equivalent to 

- fip{w) = _ 

P 

The second critical point equation dw^{0, pp{w);T, z) = 0 reduces to 

dnj{ppiw),T) = 0 . 

Since ■^pp\.^^o is always tangent to the face of P at /rp(w°), this is just the condition that t is orthogonal to 
the face of P containing p,p(w), which is automatically satisfied when r is in the normal cone to P at fj,p(w). 
Hence, in order to have critical points, r must equal Tz as before, and we obtain a critical submanifold 

Cz := {(0, w) : pp{w) = q{z)} . 

Note that dimC^ = codimF, where F is the face of P containing q{z). 

By a calculation using the methods of |SZ2[ §3.2.2] (where a more complicated phase function was used), 
one can show that the normal Hessian is nondegenerate whenever 2; is not a transition point. The asymp¬ 
totic expansion of Theorem 10 then follows by the method of stationary phase with nondegenerate critical 
submanifolds. 

Remark: We may also express polytope characters through the Szego kernel H^^ of F[^{Mp, Lp) by means 
of the obvious identity 

(176) 


XNp{Jn = / (e*‘" • w, w) dVolMp (w). 


> M p 
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In an article with T. Tate |STZ1| . we describe an explicit construction of the Szego kernel of a toric variety, 
which we use to obtain a formula for the polytope character XArp(e®‘^) as a complex oscillatory integral over 
the toric variety Mp. 
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(v) 

(Ell 

7p 

® 

7p|P 

O 

Kn 
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M 

m 


(Pii. (pii 

n|p 

(Ell. (El) 

E 

mi 

Tz 

(inii-dTKii 

Xa 

O 

IIXcll 

(El 

CP 

Aol 

(E|i. (El 

Xnp 

(El 


(E3i 


(|M]i 


(P^ 

wfs 

(El 

Ap 

® 

bp 

oni 

b{x} 

(Ell, (El) 

Cp 

(El 

Cx = CP 

6U 

E\np 

m 


)imi 


(El 

Flow(a;) 

El 

I/WIfs 

El 

In, In 

(IHHl 

J{x) 

(El 


(El 

C 

(PI 

Log 

(ITH71 

Mx 

El 

N{Q) 

jEE 

Pf 

® 

P[h) 

(EOl 

Poly(P) 

0 

q{^) 

(El^CSl 

Tip 

El 

S{e) 

(El 

Sf 

® 

rjrm 

El 

Todd(J^, d/dh) 

m 

UN 

(E^ 

Uoo 

(113511 


El 

\Zh,...,h\ 

(EOl 


(PI 
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